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SUMMARY 


Universal  curves  are  presented  to  facilitate  the  prediction  of  magnetohydro- 
dynamic  control  forces  in  rocket  motor  practice.  A  set  of  equilibrium  calculations 
also  required  for  this  purpose  is  presented  for  21  rocket  propellant  combinations 
with  alkali  atom  seeding.  Since  the  nozzle  expansion  process  can  lead  to  electri¬ 
cal  conductivities  at  the  exit  section  much  larger  than  those  corresponding  to 
local  ionization  equilibrium,  a  graphical  method  is  introduced  and  discussed  M4iich 
enables  quantitative  estimates  to  be  made  of  the  extent  of  this  nonequilibrium  en¬ 
hancement  as  well  as  the  conditions  under  which  it  can  be  ignored.  By  combining 
this  material  it  is  possible  to  extract  general  information  as  to  the  independent 
factors  favoring  maximum  exhaust  jet  interaction  with  the  applied  magnetic  field. 
Some  of  the  major  factors  (other  than  large  field  strength)  giving  rise  to  favor¬ 
able  conditions  are  a)  presence  of  low  ionization  potential  materials  in  "hot" 
chambers,  b)  small  values  of  the  effective  specific  heat  ratio,  c)  moderate  rather 
than  extremely  large  chamber  pressures,  and  d)  large  physical  size.  It  appears, 
from  the  calculations  carried  out  thus  far,  that  control  forces  of  the  order  of 
0.1  to  1.0  percent  of  the  gross  thrust  are  attainable  by  applying  large  magnetic 
fields  (10,000  to  30,000  gauss)  to  the  exhaust  jets  of  rocket  motors  in  the  10,000 
lb  thrust  class,  which  have  been  seeded  with  alkali  atoms  to  the  extent  of  0.1  to 
1.0  percent  by  mass.  With  the  use  of  superconducting  magnet  coils,  larger  field 
strengths  are  possible  (perhaps  up  to  100,000  gauss)  and  side  thrusts  of  the  order 
of  several  percent  of  the  gross  thrust  are  predicted.  Because  of  competing  fac¬ 
tors,  these  control  force  expectations  generally  apply  to  both  the  equilibriiwi  as 
well  as  nonequilibrium  cases  envisioned. 

ACKNOWLEDGEMENTS 


The  author  wishes  to  thank  Dr.  C.  Boyars  for  suggesting  the  problem  and 
Dr.  H.  F.  Calcote  for  providing  the  opportunity  to  discuss  many  of  its  controver¬ 
sial  aspects.  He  is  also  particularly  indebted  to  F.  Kuehner  for  carrying  out 
the  bulk  of  the  numerical  calculations  presented  herein,  and  to  the  many  authors 
whose  work  is  implicitly  or  explicitly  represented  in  the  bibliography  for  pro¬ 
viding  the  groundwork  iMiich  made  this  study  possible. 


ii 


TN-41 


TABLE  OF  CONTENTS 

Page 

SUMMARY .  ii 

ACKNOWLEDC3EMENTS .  ii 

LIST  OF  TABLES .  v 

LIST  OT  FIGURES .  v 

NOMENCLATURE .  vi 

I.  INTRODUCTION  AND  OBJECTIVES .  1 

II.  EFFECT  OF  AN  INCLINED  MAGNETIC  FIELD  ON  TRANSVERSE  THRUST  ....  1 

III.  THE  ELECTRICAL  CONDXTIVITY  OF  SEEDED  COMBUSTION  GASES .  2 

IV.  ANTICIPATED  VALUES  OF  THE  MAGNETOGASDYNAMIC  INTERACTION 

PARAMETER  Q .  8 

1.  Establishment  of  the  Reference  Value,  of  the  Magnetogas- 

dynamic  Interaction  Parameter  .  10 

2.  Enhancement  of  the  Interaction  Parameter  Due  to  Electron-Ion 

Recombination  Lags  in  the  Nozzle .  13 

3.  Effect  of  Physical  Size  (or  Thrust  Level)  on  the  Magnitude  of 

the  Magnetogasdynamic  Interaction  Parameter  .  17 

4.  Effect  of  Nozzle  Expansion  Conditions  on  the  Magnitude  of  the 

Magnetogasdynamic  Interaction  Parameter  at  the  Nozzle  Exit 
Section .  18 

V.  NUMERICAL  EXAMPLE;  TRANSVERSE  THRUST  FOR  A  300  PSIA  CHAMBER 
PRESSURE,  10,000  POUND  THRUST  OXYHYDROGEN  MOTOR;  0.1)6  POTASSIUM 
SEED .  19 

1.  Determination  of  Q  . .  19 

ref 

2.  Effect  of  Electron-Ion  Recombination  Lag  .  19 

3.  Effect  of  Magnetic  Field  Strength  .  19 

4.  Effect  of  Physical  Scale  (or  Thrust  Level)  .  19 

5.  Determination  of  the  Equilibrium  Expansion  Factor,  E  .  20 

VI.  DISCUSSION .  22 

1.  Effect  of  0/T^  on  Interaction  Parameter  .  23 

2.  Effect  of  Specific  Heat  Ratio .  23 

3.  Effect  of  Expansion  Pressure  Ratio  for  a  Class  of  Equal  Thrust 

Motors  Expanding  to  the  Same  Back  Pressure .  25 

4.  Scale  Effect  for  Motors  of  the  Same  Pressure  Level .  25 

5.  Accuracy  of  Physical  Model . 26 

iii 


TN-41 


TABLE  OF  COlffENTS  (CON'T) 

Page 

VII.  CONCLUDING  REMARKS .  27 

VIII.  REFERENCES .  28 

APPENDICES 

A  TRAJECTORY  ANALYSIS:  MAGNETOGASDYNAMIC  JET  DEFLECTION  .  31 

B  ELECTRON- ION  RECOMBINATION  LAGS  .  34 

Formulation .  34 

"Sudden  Freeze"  Approximation  .  36 

Calculation  Procedure  .  37 

Selection  of  Realistic  Parameters  .  39 

Accuracy  of  the  Sudden  Freeze  Approximation;  Improved  Approx¬ 
imations  . 

C  DECAY  OF  ELECTRICAL  CONDUCTIVITY  DOWNSTREAM  OF  THE  NOZZLE  ....  45 


FIGURES 


48 


TN-41 


LIST  OF  TABLES 

Table  Pa9e 

I  DATA  ON  THE  SINGLE  IONIZATION  OF  SEVERAL  ELEMENTS  AND  MOLECULES  ...  4 

II  VALUES  OF  -  LOGio  FOR  THE  SINGLE  IONIZATION  OF  ALKALI  ATOMS  ...  7 

III  ROCKET  MOTOR  DATA  FOR  CHAMBERS  OPERATING  AT  300  PSIA .  11 

IV  ROCKET  MOTOR  DATA  FOR  CHAMBERS  OPERATING  AT  1000  PSIA .  12 

V  TEMPERATURE  AND  PRESSURE  DEPENDENCE  OF  THE  ION  RECOMBINATION 

COEFFICIENT  (ABRIDGED)  .  15 

VI  EFFECTIVE  VALUES  OF  d(ln  cj)/d(ln  T)  FOR  CESIUM  SEEDING  IN  ROCKET 

MOTOR  APPLICATIONS .  24 

VII  PRESSURE  AND  TEMPERATURE  DEPENDENCE  OF  THE  ION  RECOMBINATION 

COEFFICIENT .  40 

VIII  THE  FUNCTION  H(M)  FOR  Y  =  1.2  ISENTROPIC  EXPANSIONS  IN  A  <t>  =  15° 

HALF  ANGLE  CONE .  44 

LIST  OF  FIGURES 

Figure  Page 

1  GEOMETRIC  CONFIGURATION:  IONIZED  JET  PASSING  THROUGH  AN  IMPOSED 

MAGNETIC  FIELD .  48 

2  DEPENDENCE  OF  ATTAINABLE  TRANSVERSE  THRUST  ON  THE  MAGNITUDE  OF  THE 

MAGNETOGASDYNAMIC  INTERACTION  PARAMETER  Q .  49 

3  CORRELATION  OF  Q  .  WITH  CHAMBER  TEMPERATURE  FOR  300  PSIA  MOTORS  .  .  50 

re  I 

4  THE  EFFECT  OF  ELECTRON-ION  RECOMBINATION  LAG  ON  EXIT  VALUES  OF  THE 
ELECTRICAL  CONDUCTIVITY  FOR  ISENTROPIC  NOZZLE  EXPANSIONS 

(Y  =  1.2|  Langevin  Recombination) .  51 

5  THE  EFFECT  OF  ELECTRON-ION  RECOMBINATION  LAG  ON  EXIT  VALUES  OF  THE 
ELECTRICAL  CONDUCTIVITY  FOR  ISENTROPIC  NOZZLE  EXPANSIONS 

(Y  =  1<2(  Bates  Recombination) .  52 

6  THE  EFFECT  OF  ELECTRON-ION  RECOMBINATION  LAG  ON  EXIT  VALUES  OF  THE 
ELECTRICAL  CONDXTIVITY  FOR  ISENTROPIC  NOZZLE  EXPANSIONS 

(Y  =  1.2j  Thompson  Recombination) .  53 

7  NOZZLE  SHAPE  TREATED  IN  RECOMBINATION  STUDY  .  54 

8  GRAPHICAL  DETERMINATION  OF  THE  EFFECT  OF  IONIZATION  NONEQUILIBRIUM 
ON  EXIT  VALUES  OF  THE  ELECTRICAL  CONDUCTIVITY  FOR  ISENTROPIC 

NOZZLE  EXPANSIONS  .  55 

9  PRESSURE  DEPENDENCE  OT  EQUILIBRIUM  EXPANSION  FACTOR*  (Y  =  1.2)  ...  56 

10  PRESSURE  DEPENDENCE  OF  EQUILIBRIUM  EXPANSION  FACTOR;  (Y  =  1.3)  ...  57 

11  PRESSURE  DEPENDENCE  OF  THE  ION  RECOMBINATION  COEFFICIENT  FOR  AIR  .  .  58 


V 


TN-41 


LIST  OF  FIGURES  (CON'T) 


Figure  Page 

12  LCKATION  OF  THE  MACH  NUMBER  BEYOND  WHICH  THE  ELECTRON  MOLE 
FRACTION  IS  EXPECTED  TO  REMAIN  UNCHANCSD 

(Y  =  1.2j  Thompson  Recombination) .  59 

13  EXACT  INTEGRAL  CURVES  OF  THE  DIFFERENTIAL  EQUATION  GOVERNING 

CONDUCTIVITY  CHANGES  IN  THE  NOZZLE .  60 

14  IMPROVED  APPROXIMATE  SOLUTIONS;  MAXIMUM  LOWER  BOUND  TECHNIQUE  ...  6l 

15  COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTIONS  FOR  A  PARTICULAR  CASE 

(0Aj.  =  15,  =  10-^,  Y  =  1.2,  =  -  7/2,  u)p  =  +  1) .  62 

16  MATCH  POINTS  ON  THE  LOG  5  -  LOG  p-’  PLANE;  MAXIMUM  LOWER  BOUND 

METHOD  (y  =  1.2) .  63 

17  DECAY  OF  ELECTRICAL  CONDUCTIVITY  DOWNSTREAM  OF  THE  NOZZLE  .  64 

18  DECAY  OF  SPATIALLY  AVERAGED  ELECTRICAL  CONDUCTIVITY  DOWNSTREAM 

OF  THE  NOZZLE .  65 

19  UNIVERSAL  FUNCTION  REQUIRED  TO  CALCULATE  THE  STREAMWISE  VARIABLE 

z  for  GIVEN  VALUES  OF  C  and  x/d  .  66 

r  ex 


NOMENCLATURE 

A  cross-sectional  area  at  a  station  in  the  nozzle 

I  A/A^  =  (d/d^)2 

B  magnetic  field  strength 

c  ideal  exhaust  velocity 

c’^  characteristic  velocity  of  propellant 

C  ion  recombination  rate  parameter;  Eq.  (l8) 
r 

Cs  chemical  symbol  for  the  element  cesium 

d  diameter  of  cross-section 

D^  throat  diameter  (=  d^) 

e  electronic  charge 

energy  change  across  ionization  reaction 
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h  Planck's  constant 
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k  Boltzmann's  constant 


vi 


TN-41 


K 

P 

K 

L 

M 

m 

ifi 

% 

n 

n 

P 

P 

Q 

R 

S 

T 

f 

u 

u 

'fi 

X 

X. 

1 

X 

Y 

y 


Greek 


a 

a 

] 

Y 

e 

e 

p 

a 


NOMENCLATURE  (CON’T) 

equilibrium  constant}  Eq,  (6) 
chemical  symbol  for  potassium 

effective  length  of  the  magnetic  field  in  the  axial  direction 

Mach  number 

mass  of  the  electron 

mass  flow  rate 

molecular  weight 

number  density 

n/nc 

total  pressure 

p/Pc 

magnetogasdynamic  interaction  parameter;  Eq,  (l) 

effective  collision  cross-section 
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represents  seed  species 

absolute  temperature 

streamwise  velocity 
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mass  rate  of  production  per  unit  volume  (source  term)}  Eq.  (15) 

nondimensional  axial  position  [Eq.  (5)] 

mole  fraction  of  species  i  in  mixture  [Eq.  (2)] 

distance  in  axial  direction 

mass  fraction  [Eq.  (B-6)]  or  nondimensional  transverse  position  [Eq.  (A-5) 
distance  in  transverse  direction 


mass  fraction  of  seed  material  added 
ion  recombination  coefficient}  Eq.  (l5) 
specific  heat  ratio  of  carrier  gas  (y  = 

inclination  angle  between  imposed  magnetic  field  and  nozzle  axis 
characteristic  temperature  for  ionization}  Table  I 
total  density 

electrical  conductivity}  Eq.  (2) 
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NOMENCLATURE  (CON’T) 

£  nonequilibrium  enhancement,  Eq.  (C-3) 

T  characteristic  time  or  nondimensional  timej  Eq.  (A-5) 
divergence  half  angle  of  conical  expansion  section 
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angle  parameter  in  Eqs.  (B-15),  (B-16) 

Subscripts 

+  pertaining  to  positive  ions 

X  in  transverse  direction 

0  undisturbed  (without  magnetic  field);  at  time  =0;  or  at  some  arbitrary 
initial  point  in  nozzle  [see  Eq.  (B-27)] 

*  pertaining  to  sonic  conditions  (Mach  number  =  unity) 

c  in  the  combustion  chamber 
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I.  INTRODUCTION  AND  OBJECTIVES 

An  electrically  conducting  gas  can  be  subjected  to  a  force  when  it  moves 
through  a  stationary  magnetic  field  since  electrical  currents  are  induced  which 
can  interact  with  the  imposed  field.  This  being  the  case,  it  is  possible,  at  least 
in  principle,  to  use  magnetic  fields  to  control  the  motion  of  high  speed  jets,  as, 
for  example,  would  be  required  in  thrust  vectoring  a  propulsion  device.  Assuming 
such  a  technique  would  have  advantages  over  more  conventional  control  methods,^  it 
becomes  necessary  to  inquire  into  the  nature  of  the  interaction  and  particularly 
the  magnitude  of  the  magnetic  field  strengths  required  to  achieve  the  desired  con* 
trol  forces.^  The  results  of  an  investigation  of  this  type  are  included  in  the 
present  work,  where  emphasis  will  be  given  to  their  presentation  in  a  form  having 
sufficient  generality  to  be  adaptable  to  different  design  needs.  Thus,  insofar  as 
possible,  the  problem  will  be  divided  into  independent  aspects  governing  the  mag¬ 
nitude  of  the  anticipated  MHD  effect  for  the  case  of  conventional  rocket  motors  ex¬ 
panding  the  products  of  combustion  reactions.  Calculations  for  a  number  of  repre¬ 
sentative  propellant  combinations  will  be  included. 

II.  EFFECT  OF  AN  INCLINED  MAGNETIC  FIELD  ON  TRANSVERSE  THRUST 

We  first  imagine  a  unifoxm  jet  of  electrically  conducting  gas  passing 
through  an  imposed  field  of  strength  B  inclined  at  an  angle  6  to  the  direction  of 
flow  (see  Fig.  1).  If  the  jet  has  a  diameter  d  and  the  magnetic  field  has  an  ax¬ 
ial  extent  L,  then  a  fluid  particle  trajectory  analysis  (see  Appendix  A)  can  be 
carried  out  to  determine  the  anticipated  transverse  thrust  obtained  for  pre¬ 
scribed  values  of  the  field  length  to  jet  diameter  ratio,  L/d,  and  the  dimension¬ 
less  Biagnetogasdynamic  interaction  parameter 

(1) 


^Important  considerations  are  mechanical  complexity,  weight,  response,  severity 
of  high  temperature  (materials)  problem. 

^Recent  developments  in  solid  state  physics  suggest  that  superconducting  magnet 
coils’  may  become  available  which  would  enable  large  increases  in  field  strengths 
over  those  currently  attainable. 
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Here  d  is  the  electrical  conductivity  of  the  partially  ionized  gas.  The  product 
of  the  gas  density  p  with  the  exhaust  velocity  c  appearing  in  the  denominator 
will  be  recognized  as  momentum  of  the  gas  per  unit  volume.  If  F  designates  the 
streamwise  thrust  in  the  absence  of  magnetic  field,  then  for  a  fixed  field  in¬ 
clination  0  the  value  of  indeed  depends  only  upon  L/d  and  Q.  This  relation 

is  shown  in  Fig.  2  for  the  field  inclination  (0  =  45°)  giving  maximum  transverse 
thrust.  Since  the  value  of  L/d  is  expected  to  be  determined  principally  by  geo¬ 
metric  and  mechanical  considerations,  it  is  seen  that  the  determination  of 
reduces  to  the  determination  of  realistic  values  for  the  interaction  parameter  Q 
for  rocket  nozzle  exhaust  Jets.  This  is  described  in  succeeding  sections. 

III.  THE  ELECTRICAL  CONDUCTIVITY  OF  SEEDED  COMBUSTION  GASES 

The  electrical  conductivity,  d,  which  enters  the  expression  for  the  Inter¬ 
action  parameter  Q  is  consequence  of  the  presence  of  mobile  charge  carriers  in 
the  gas;  i.e.,  the  gas  must  be  partially  ionized  to  become  a  moderately  good  con¬ 
ductor.  When  the  effects  of  electron-neutral  (short  range)  collisions  dominate 
those  of  the  coulomb  interactions  (long  range)  between  electrons  and  positive 
ions,  the  electrical  conductivity  of  a  partially  ionized  gas  is  simply  related 
to  the  prevalent  mole  fraction  X  of  electrons.  Kinetic  theory  provides  the  use¬ 
ful  result* 


a  =  0.532 


(m^  kT)^ 


(2) 


where  is  the  effective  electron-neutral  collision  cross-section,  e  and  m  are 
o  e 

the  charge  and  mass  of  the  electron,  respectively,  and  T  is  the  gas  temperature. 

In  MKS  units,  if  we  introduce  the  numerical  values 

e  =  1.6  X  10" ’’  coulomb 
m^  =  0.911  X  10-*®  kg 
k  =  1.38  X  10"**  kg(meter)* (sec)"* (®K)“ ‘ 

and  express  temperature  in  degrees  Kelvin,  then  the  units  of  electrical  conduc¬ 
tivity  are  mhos/meter.  Thus  the  electrical  conductivity  can  be  estimated  if  the 


^It  will  be  assumed  that  the  gas  temperature  and  electron  temperature  are  every¬ 
where  equal. 
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gas  temperature*  T*  electron  mole  fraction  X  ,  and  collision  cross-section 
are  known.  For  the  purposes  of  the  present  study,  it  will  be  assumed  that  the 
value  of  the  electron-neutral  cross-section,  Q^,  will  not  be  very  different  from 
10"’ 9  (meter)^  for  most  combustion  gases^  of  interest  and  that  this  cross-section 
has  no  appreciable  dependence  on  temperature  level^  (see  Section  VI,  Part  5). 

It  will  further  be  assumed  that  the  condition  of  thermochemical  equili¬ 
brium  exists  for  the  combustion  gases  at  the  nozzle  inlet  section.  This  means 
that  the  electron  mole  fraction  corresponding  to  chamber  conditions  is  de¬ 
termined  by  thermodynamic  considerations  alone.  Since  the  electrons  present  in 
the  gas  must  originate  from  ionization  of  some  of  the  molecules  present  in  the 
chamber,  it  is  clear  that  X^  will  be  large  if  an  easily  ionizable  species  is 
present  in  the  product  gases. 

Table  I  shows  values  of  the  first  ionization  potential*  for  a  selected 
group  of  atoms  and  molecules  of  interest  in  the  present  discussion.  Also  tabu¬ 
lated  are  the  corresponding  values  of  the  energy  change  A  E° across  the  ioniza¬ 
tion  reaction  (in  kcal/mole)  and  the  corresponding  "characteristic  ionization 
temperature"  0  =  AE°/R,  where  R  is  the  universal  gas  constant  (1.9872  cal/deg- 
mole).  This  temperature  parameter  is  of  physical  interest  since  equilibrium  sta¬ 
tistical  mechanics  reveals  that  the  extent  of  ionization  will  necessarily  be 
small  if  the  absolute  gas  temperature  is  small  compared  with  0.  Now,  in  rocket 
motor  practice  the  chamber  temperature  levels  are  in  the  range  of,  say,  2000°K  to 
4000°K.  But  for  the  molecules  which  would  be  expected  to  prevail  in  ordinary 
combustion  products  (viz.  H^O,  O3,  CO2,  CO,  etc.)  values  of  the  characteristic 
temperature  for  ionization  are  in  the  range  of  145,000°K  to  l67,000®K.  By  com¬ 
parison,  however,  elements  exist  for  which  0  1?  .'ubstantially  lower.  The  class 
of  elements  having  the  lowest  values  of  the  characteristic  ionization  temperature 
0,  is  the  alkali  group,  with  cesium  displaying  r.  ©  of  only  44,900‘’k. 

It  can  be  verified  a  posteriori  that  the  conclusion  one  is  tempted  to  draw 
from  this  comparison  is  indeed  valid,  i.e. ,  if  one  does  not  intentionally  "seed” 
the  usual  products  of  combustion^  with  an  easily  ionizable  substance,  such  as 


^AE^  (in  kcal/mole)  can  be  obtained  from  tabulated  values  of  the  ionization  poten 
tial  (in  electron  volts)  by  multiplication  with  23.0631. 

^This  may  not  be  true  for  solid  propellants  containing  metallic  salts  or  metallic 
additives,  such  as  aluminum. 
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TABLE  I 

DATA  ON  THE  SINGLE  IONIZATION  OF  SEVERAL  ELEMENTS  AND  MOLECULES 


Species 

% 

Ionization  Energy  Parameter 

V,ev. 

AE°,Kcal/mole 

0,°K 

Lithium 

6.940 

5.363 

123.7 

62,240 

Sodium 

22.991 

5.12 

118. 

59,400 

Potassium 

39.100 

4.318 

99.59 

50,110 

Rubidium 

85.48 

4.159 

95.92 

48,270 

Cesium 

132.91 

3.87 

89.3 

44,900 

Aluminum 

26.98 

5.96 

137. 

69,200 

Calcium 

40.08 

6.09 

140. 

70,700 

Water  Vapor 

18.016 

12.56 

289.7 

145,800 

Oxygen  (O2) 

32.000 

12.5 

288. 

145,000 

Carbon  Dioxide 

44.011 

14.4 

332. 

167,000 

Carbon  Monoxide 

28.011 

14.1 

325. 

164,000 

Ammonia 

17.032 

11.2 

258. 

130,000 

Nitric  Oxide 

30.008 

9.5 

220. 

110,000 
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potassium  or  cesium  (or  a  compound  thereof)  then  the  electron  mole  fraction,  and 
hence  the  electrical  conductivity  of  the  gas  mixture,  will  be  inconsequential 
(insofar  as  the  present  magnetogasdynamic  application  is  concerned).  Conversely, 
if  one  does  "seed"  with  such  a  material  the  electrons  which  one  will  find  in  the 
gas  will, overwhelmingly, be  offspring  of  the  seed  material.  In  this  simple  case 
one  can  then  calculate  the  electron  mole  fraction  by  considering  a  single  re¬ 
versible  reaction,  written  symbolically  as 

S  Z  +  e  (3) 

Now,  in  practice,  only  small  amounts  of  the  seed  material  S  would  be  added  to  the 
combustion  mixture  to  achieve  the  desired  result.’^,®  If  S  is  added  with  a  mass^ 
fraction,  a^,  then  the  electron  mole  fraction  is  found  from  the  law  of  mass 
action  to  be  proportional  to  the  square  root  of  o^,  i.e. 

^  =  (Kp/P)^  (a^)^  (4) 

where  Kp(T)  is  the  equilibrium  constant  for  the  reaction  (3)  and  the  script  Ul’s 
represent  molecular  weights.  Combining  Eqs.  (2)  and  (4)  with  the  assumption 
=  constant,  we  find  the  functional  dependence 

a  OC  p"^  CKp(T)]^  (5) 

In  this  relation,  the  temperature  (dependence  of  the  equilibrium  constant  can  be 
obtained  from  the  Saha  equation  of  equilibrium  statistical  mechanics.^* For 
the  case  of  single  ionization  of  alkali  atoms^^  the  appropriate  relation  is 


This  is  much  more  convenient  than  the  more  commonly  used  mole  fraction  for  spec¬ 
ifying  the  seed  concentration. 

^subject  to  the  assumption  that  the  seed  material  does  not  undergo  appreciable 
reaction  with  the  combustion  products 

this  case  the  pre-exponential  factor  giving  the  dependence  on  electronic  par¬ 
tition  function  ratios  becomes  approximately  1.00  in  the  temperature  range  of 
interest  here.  This  factor  has  therefore  been  omitted  from  Eq.(6). 
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3  5 

Kp(T)  =  i2mn^)^  h-3  (kT)^  exp  (-  ©A)  (6) 

where  h  is  Planck’s  constant  and  all  other  quantities  have  already  been  defined. 
Therefore,  subject  to  these  approximations,  the  electrical  conductivity,  cJ,  of 
a  seeded  mixture  of  combustion  gases  will  depend  upon  pressure  and  temperature 
in  accord  with  the  law” 


d 


exp  (-^©A) 


(7) 


This  relation  will  be  used  to  describe  the  variation  in  electrical  conductivity 
for  equilibrium  isentropic  nozzle  expansions  of  seeded  combustion  gases. 

Values  of  the  equilibrium  constant  computed  from  Eq.  (6)  are  given  in 

Table  II  for  the  useful  range  of  temperatures  and  ©A«  Using  Table  II,  one  can 

rapidly  evaluate  the  equilibrium  constant  Kp  for  any  chamber  temperature  and 

seed  material  (i .e. ,  any  ©).  The  electron  mole  fraction  is  then  evaluated  from 

Eq.(4)  and  inserted  into  the  electrical  conductivity  relation  (2).  It  is  seen 

that  the  absolute  value  of  conductivity  at  the  nozzle  inlet  section  is  therefore 

determinable  from  a  knowledge  of  the  gas  pressure,  p  ,  temperature,  T  ,  and  mean 

c  c 

molecular  weight,  Hl^,  together  with  the  nature  of  the  seed  material  (©,  111^)  and 
the  fractional  extent  of  seeding  (Og).  If  equilibrium  were  maintained  during  the 
nozzle  expansion  the  electrical  conductivity  would  vary  in  accord  with  Eq.  (7), 
since  the  mole  fraction  of  ionizable  material  would  be  sensibly  constant. 
Moreover,  the  seed  material  may  often  be  regarded  as  being  immersed  in  a  carrier 
gas  which  expands  isentropically  through  the  nozzle.  Then,  if  the  chemical  com¬ 
position  of  this  carrier  gas  does  not  change,^  pressure  and  temperature  in  the 
nozzle  will  be  related  in  accord  with” 


Y 


(8) 


It  follows  that  if  ionization  equilibrium  were  maintained  everywhere,  then  for  a 
prescribed  value  of  the  specific  heat  ratio,  y,  the  conductivity  ratio  o/o 

c 


t 


i.e.,  if  the  expansion  process  is  chemically  frozen 
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TABLE  II 

VALUES  OF  -  LOG,o  FOR  THE  SINGLE  IONIZATION  OF  ALKALI  ATOMS 


T 

eAc 

T  K 

0 

10 

20 

30 

40 

1000 

-1.0172 

3.3257 

7.6686 

12.0116 

16.3545 

1500 

-1.4575 

2.8855 

7.2284 

11.5714 

15.9143 

2000 

-1.7698 

2.5731 

6.9161 

11.2590 

15.6020 

2500 

-2.0121 

2.3308 

6.6738 

11.0167 

15.3597 

3000 

-2.2101 

2.1329 

6.4758 

10.8188 

15.1617 

3500 

-2.3774 

1.9655 

6.3085 

10.6514 

14.9944 

4000 

-2.5224 

1.8205 

6.1635 

10.5064 

14.8494 

4500 

-2.6503 

1.6927 

6.0356 

10.3786 

14.7215 

5000 

-2.7647 

1.5783 

5.9212 

10.2642 

14.6071 
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would  be  a  function  only  of  the  pressure  ratio  p  /p  (say)  and  the  parameter 

c 

©Aj.«  In  connection  with  a  subsequent  nozzle  flow  recombination  study*  this 

conductivity  ratio  function  will  be  shown  graphically  in  Figs,  4*  5  and  6»  for 

a  carrier  gas  specific  heat  ratio  of  1.2.  The  values  of  o/o  can  now  be  brought 

c 

to  bear  on  the  determination  of  the  magnetogasdynamic  interaction  parameter  Q 
for  the  exhaust  section  of  rocket  motors, 

IV.  AhTTICIPATED  VALUES  OF  THE  MAGNETOGASDYNAMIC  INTERACTION  PARAMETER  Q 

While  we  are  interested  in  the  value  of  the  magnetogasdynamic  interaction 
parameter  Q  at  the  nozzle  exit  section,  a  cursory  study  of  the  rocket  motor 
literature' reveals  that  most  information  for  various  propellant  combina¬ 
tions  is  tabulated  at  chamber  conditions.  Thus  for  prescribed  chamber  pressures 
(usually  300  psia*  or  1000  psia)  extensive  data  is  available  giving  mean  molecular 
weight*  chamber  temperature,  ratio  of  specific  heats*  characteristic  velocity  (c*) 
and  sometimes  detailed  species  concentrations.  This  abundance  of  chamber  data 
suggests  that  for  design  purposes  it  will  be  convenient  to  split  the  determination 
of  the  interaction  parameter  Q  at  the  nozzle  exit  section  into  several  more  or 
less  independent  factors: 

(a)  dependence  on  particular  propellant/additive  system 

(b)  dependence  on  chemical  kinetic  parameters  governing  departures 
from  ionization  equilibrium^ 

(c)  dependence  on  available  magnetic  field  strength 

(d)  dependence  on  physical  scale  (thrust  level)  other  than  that  due  to  (b) 

(e)  dependence  on  expansion  conditions  (pressure  ratio  across  nozzle) 

This  can  be  accomplished  by  rewriting  Eq.  (l)  as  follows 


Throughout  this  work  the  term  nonequi librium  deals  with  the  inability  of  the 
electron-ion  recombination  rate  to  maintain  equilibrium  in  an  expanding  flow. 
This  should  be  distinguished  from  the  use  of  the  phrase  nonequilibrium  ioniza¬ 
tion  to  describe  such  phenomena  as  chemi-ionization.*^ 
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Q 


=  Q 


ref 


(9) 


where 


j  =  ^ref  ‘^ref 

^ref  “  pc* 


(10) 


and  c*  is  the  characteristic  velocity  of  the  propellant.  The  factors  in  this 
decomposition  have  the  following  significance; 


(a)  fixes  the  level  of  the  magnetogasdynamic  parameter  for  some  ref¬ 
erence  motor  size  and  field  strength  and  is  computed  based  on  chamber 

conditions.  As  will  be  seen,  Qj>gf  can  be  calculated  for  each  propellant  system 
for  which  theoretical  performance,  chamber  temperature,  and  molecular  weight  data 
are  available  and  will  depend  upon  the  nature  and  amount  of  the  "seed*’  material 
added  to  enhance  the  electrical  conductivity.  In  the  MKS  system  of  units  con¬ 
venient  choices  for  and  d^^^  might  be  1  weber/meter*  and  1  meter,  respec¬ 
tively.  Numerical  values  for  in  practical  cases  will  be  given  later  on. 


(b)  The  factor  d  /p  _  represents  the  enhancement  in  the  conductivity 

©X  ©X  9  ©CJ 

at  the  exhaust  section  which  would  result  from  the  inability  of  ion-electron  re¬ 
combination  rates  to  keep  pace  with  the  rapidly  changing  pressure-temperature 
environment  in  the  nozzle. This  enhancement  will  be  discussed  in  some  de¬ 
tail  in  Part  2  of  this  section,  as  well  as  Appendix  B. 


(c)  The  factor  (b/B^^^)^  accounts  for  the  effect  on  the  interaction 
parameter  Q  of  field  strength  capabilities  differing  from  Bj.gf 

(d)  The  factor  d./d  ,,  being  the  ratio  of  the  actual  throat  diameter 

t,  r©t 

to  d  .a  accounts  for  the  effect  on  Q  of  the  scale  of  the  motor  In  question.  As 
ref  ■ 

will  be  seen  in  Part  3  of  this  section  d./d  -  is  proportional  to  the  square 

t  ret 

root  of  the  thrust  level  for  a  given  chamber  pressure.  Scale  effects  are  also 
implicit  in  part  (b)  above,  but  these  require  separate  treatment. 


(e) 


In  Eq. 


(9)  the  factor  in  brackets,  i .e. 

\ 


P  u 


(11) 
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is  a  function  of  only  3  variables:  0/T  >  specific  heat  ratio,  y*  expansion 

c 

pressure  ratio,  p^/p.  It  is  precisely  that  function  by  which  Qpgf  (°'ex^'^ex,eq^ 

(d^/d^g^)  must  be  multiplied  to  recover  Q  =  ^/(pgj^  c).  Discussion 

of  this  "equilibrium  expansion  factor",  E,  will  be  postponed  to  Part  4  of  this 
section. 

On  the  basis  of  this  decomposition  it  becomes  possible  to  combine  calcu¬ 
lations  of  for  all  of  the  usual  propellant  combinations  with  data  read  from 

a  few  universal  graphs  to  arrive  at  the  value  of  the  interaction  parameter  Q,  ex¬ 
pected  in  the  jet  exhaust  for  any  desired  expansion  ratio,  seed  material  and  con¬ 
centration,  magnetic  field  strength  and  physical  scale.  A  knowledge  of  Q  and  the 
effective  field  length  (measured  in  jet  diameters)  is  then  all  that  is  needed  to 
determine  the  control  force  (Fig.  2). 

1.  Establishment  of  the  Reference  Value,  of  the  Magnetogasdynamic  Inter- 

action  Parameter 


Since  it  will  be  seen  that  it  is  easy  to  correct  the  results  of  calcula¬ 
tions  of  Q  ,  for  the  effects  of  changing  either  the  seed  material  or  the  seed 
ref 

mass  fraction,  numerical  results  will  be  given  here  for  a  class  of  rocket  motors 
operating  with  chamber  pressures  of  300  psia  and  1000  psia,  respectively,  seeded 
with  cesium  to  the  extent  of  one  percent  by  mass.  For  this  case:  0  =  44»900*^K, 

=  10"*  and  111^  =  132.91.  Reference  quantities  are  =  1  meter  and  B^^^  =  1 

Weber/meter*  =  10,000  gauss.  Numerical  results  for  seven  fuel/oxidizer  combina¬ 
tions  at  300  psia  chamber  pressure’^*’*’’’  are  collected  in  Table  III.  Results 
for  14  fuel/oxidizer  combinations  at  1000  psia  chamber  pressure’^**®  are  collec¬ 
ted  in  Table  IV. 

In  considering  the  numerical  values  of  Q^gf  tabulated  here,  it  is  found 
that  they  correlate  almost  solely  with  chamber  temperature.  This  indicates  that 
differences  in  Kl  and  y  from  propellant  to  propellant  are  dominated  by  the  expon¬ 
ential  dependence  of  electrical  conductivity  on  chamber  temperature  (see  Eq.  (7)). 
This  correlation  is  shown  in  Fig.  3,  where  is  plotted  vs.  chamber  tempera¬ 

ture  (on  log-log  paper)  and  the  individual  points  are  identified  with  the  propel¬ 
lants  appearing  in  Table  III.  This  means  that  to  a  good  approximation 

3 

Q„f  oc  T^*  exp  (-  i0/T^)  (12) 
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TABLE  III 

RXKET  MOTOR  DATA^  FOR  CHAMBERS*  OPERATING  AT  300  PSIA 


Propellants 

O/F  Ratio 

^ref 

T 

c 

c* 

c 

Yc 

Lox-Ammonia^^ 

1.4 

1.5 

X 

10-* 

3010 

1770 

2500 

1.23 

Lox-Gasolene 

2.5 

3.5 

X 

o 

1 

M 

3295 

1690 

2370 

1.22 

Lox -Hydrazine 

0.5 

0.71X 

10-2 

2750 

1845 

2542 

1.23 

Lox-Kerosene 

2.2 

4.0 

X 

10-2 

3350 

1737 

2428 

1.24 

Flourine-Hydrazine 

1.9 

30. 

X 

10-2 

U30 

2106 

2928 

1.33 

Nj04-Hydrazine 

1.0 

1.5 

X 

1— • 
o 

1 

2975 

1734 

2420 

1.27 

Hydrogen 

7.937 

4.7 

X 

o 

1 

3437 

2077 

2920 

1.20 

TABLE  III  (continued) 


Propellants 

wm 

”e,c 

Reference 

Lox-Ammonia^^ 

19.7 

0.63 

X 

10"“^ 

4.4  X 

10' 

3.1 

X 

10' 3 

13 

Lox-Gasolene 

22.7 

1.49 

X 

10-“^ 

1.0  X 

103 

0.68 

X 

10'* 

13 

Lox-Hydrazine 

15.9 

2.58 

X 

10"  5 

1.9  X 

10' 

1.4 

X 

10' 3 

13 

Lox-Kerosene 

21.6 

1.67 

X 

lO"-^ 

1.1  X 

10* 

0.74 

X 

10'* 

13 

Flourine-Hydrazine 

18.0 

1.10 

X 

10-3 

0.63X 

103 

3.7 

X 

10'* 

13 

NjO^-Hydrazine 

18.8 

0.58 

X 

10"^ 

4.1  X 

10' 

2.9 

X 

10' 3 

13 

Hydrogen-Oxygen 

15.7 

1.70 

X 

10"'^ 

1.1  X 

10* 

0.74 

X 

10'* 

19 

Units  are  as  follows:  T  in  c*  in  meter  (sec)“';  Y  dimensionless;  HI  in 

c  c  c 

gn  (gn-mole)"’j  X  dimensionless;  a  in  mho  (meter)"' ;  n  in  electrons  (cm)"^, 
c  in  meter  (sec)"' 

^  Seeded  with  1  percent  cesium  (by  mass);  ^  meter;  =  1  Weber  (meter) 

^^Lox  =  liquid  oxygen 
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and  new  values  of  (*t  the  same  chamber  pressure)  can  be  estimated  simply 

from  a  knowledge  of  the  chamber  temperature  alone»  and  any  value  of  Qj.g£  already 
given  here. 

It  is  also  relatively  simple  to  take  the  values  of  Qj,^f  computed  for  ces¬ 
ium  seeding  and  convert  them  to  apply  to  other  seeding  materials,  since  we  have 

already  seen  that  if  the  seed  mass  fraction,  a  ,  is  preserved 

s 

Qref  (13) 


Thus,  if  we  have  two  seed  materials,  designated  by  the  subscripts  1  and  2,  the 
value  of  will  change  in  accord  with 


Q 

Q 


ref  ,1 
ref  ,2 


(14) 


where  A  Q  =  0i  -  ©2. 

Corrections  to  account  for  smaller  (and  perhaps  more  realistic"^)  seed  mass 
fractions  are  likewise  made  with  the  use  of  the  square  root  relation 

oc  (Og)^.  Thus,  values  for  a  tenth  of  a  percent  (by  mass)  cesium  seed  would 
be  smaller  than  the  values  given  in  Tables  III  and  IV  by  l/'\/lO  =  .3162.  In  prac 
tice,  an  optimum  seed  concentration  will  probably  exist,  as  has  been  described  by 
Rosa‘S  and  Sherman.* 

2.  Enhancement  of  the  Interaction  Parameter  Due  to  Electron-Ion  Recombination 
Lags  in  the  Nozzle 

The  nonequilibrium  effect  is  intrinsically  the  most  difficult  to  assess 
with  any  generalityj  however,  some  progress  can  be  made  by  drawing  on  the  work 
of  Bray,*’  Hall,**  et.al.  on  deviations  from  dissociation  equilibrium  in  hyper¬ 
sonic  nozzle  expansions,  as  well  as  that  of  Smith,’*  and  Eschenroeder  and 
Daiber.’*^  In  principle,  for  a  nozzle  of  a  prescribed  shape,  a  weak  coupling 
(between  the  ion  species  and  the  carrier  gas)  approximation  allows  the  history 
of  the  reaction  to  be  readily  obtained  from  a  consideration  of  the  conservation 
equation  for  electrons  alone.  This  equation  will  contain  a  convective  teim  (see 
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Eq.  B-6)  and  a  recombination  (chemical  kinetic)  teim,^  The  latter  is  expected 
to  take  the  nonlinear  form  (see  Eq.  B-5) 

(n^,,^)»]  (15) 

where  is  the  net  mass  rate  of  production  of  electrons  per  unit  volume.  Here 
“j.  is,  formally,  the  ion  recombination  coefficient.  Actually,  the  dominant  elec¬ 
tron  loss  mechanism  need  not  be  "bimolecular"  so  that,  in  general,  the  ion  recom¬ 
bination  coefficient  appearing  in  Eq.  (15)  can  have  a  dependence  upon  total  pres¬ 
sure  as  well  as  temperature.  This  is  usually  expressed  by  a  power  law  of  the  form 

UJ  U) 

oc  T  *  p  P  (16) 

where  the  appropriate  set  of  exponents  (uj^,u)p)  depends  upon  the  anticipated  micro¬ 
scopic  mechanism.  As  described  in  Appendix  B,  three  possible  sets  of  values  are 
given  in  Table  V.  Because  of  the  extensive  variation  of  temperature,  and  particu¬ 
larly  pressure  in  a  nozzle  flow,  predictions  of  electron  recombination  lags  using 
each  of  the  models  above  differ  considerably,  as  will  be  seen. 

Our  approach  to  this  problem  has  been  to  consider  a  family  of  realistic 
nozzle  shapes  and  to  write  a  nondimensional  differential  equation  satisfied  by 
the  electrical  conductivity  itself.*  To  obtain  a  body  of  approximate  results  for 
immediate  use,  we  have  invoked  the  observations  of  Bray  and  Hall,  etal. ,  that  the 
mole  fraction  of  "relaxing"  species  is  likely  to  suddenly  "freeze"  out  to  an 
asymptotic  value  somewhere  in  the  nozzle,  it  is  possible  to  make  a  rough  estimate 
of  the  conductivity,  of  the  gas  at  the  freeze  point.  In  accord  with  Eq.  (2), 
the  electrical  conductivity,  to  be  expected  at  the  nozzle  exit  will  then  be 

given  by 

“ex  <”) 


^axial  diffusion  is  neglected 

^For  radar  attenuation  and  related  problems,  the  electron  number  density  is  of 
more  direct  interest.  This  information  can  easily  be  recovered  from  the  results 

given  in  this  report  by  everywhere  making  the  replacement  dVr  and  noting 

that  X  =  n  f/p.  Therefore  n  =  p  d/V?  . 
e  e  e 
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TABLE  V 

TEMPERATURE  AND  PRESSURE  DEPENDENCE 
OF  THE  ION  RECOMBINATION  COEFFICIENT^^*^'^*^^ 
(ABRIDGED)^ 


“t 

U) 

P 

Mechanism 

7 

■  2 

1 

Three  Body,  Thompson  (1924) 

4 

-1 

Three  Body,  Langevin  (1903) 

.  2 

2 

0 

Dissociative,  Bates  (1950) 

^See  Table  VII  for  more 

detail. 
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provided  the  temperature  dependence  of  the  cross-section  is  neglected. 

Within  a  given  class  of  nozzles,  if  the  values  of  the  following  parameters 

are  prescribed:  specific  heat  ratio  Y  of  the  carrier  gas,  the  exponents  and 

w  ,  and  0/T  }  then  exit  value  of  the  nondimensional  electrical  conductivity 
P  c 

will  depend  only  on  the  magnitude  of  a  nondimensional  recombination  para- 

meter 


C 

r 


a 

r,c 


n  d 
e,c  t 


c« 


(18) 


This  parameter  can  be  interpreted  physically  as  the  ratio  of  the  characteristic 

flow  time  d./c*  to  the  electron  half  life  (a  n  )"’  in  the  chamber.  A  dis- 
t  r»ce>c 

cussion  of  the  magnitudes  of  likely  to  be  obtained  in  practice  will  be  post¬ 
poned  to  Section  V.  It  can  be  verified  from  the  previous  relations  that  the  re¬ 
combination  parameter  will  scale  with  pressure  and  temperature  approximately  as 


u  +J-  w 

C  oc  p  P  T  ^  exp  (-  -J-GA)  (19) 

r 


for  a  fixed  fractional  seed  concentration. 

We  present  here  the  results  of  the  approximate  "sudden  freeze"  analysis 
given  in  detail  in  Appendix  B.  The  accuracy  of  this  method  is  also  discussed  in 
Appendix  B,  by  comparing  its  results  with  those  of  an  exact  solution.  The  results 
given  here  are  presented  as  contours  of  constant  recombination  parameter  on  the 
map  of  conductivity  ratio  vs.  pressure  ratio  for  isentropic  expansions  (see  Figs. 
4,3,6).  The  carrier  gas  is  assumed  to  have  a  specific  heat  ratio  of  Y  =  1.2  in 
each  case.  Thus  the  only  difference  between  Figs.  4,  5  and  6  is  in  the  choice  of 
the  exponents  and  describing  the  dependence  of  the  ion  recombination  co¬ 
efficient  on  temperature  and  pressure.  Fig.  6  corresponds  to  one  extreme 
choice  =  -  7/2,  =  +  1  (Three-body, Thompson)],  whereas  Fig.  4  corresponds  to 

another  extreme  choice  [viz,  u)^  =  +  ^,  =  -  1  (Three-body, Langevin)].  Fig.  5, 

corresponding  to  =  -  3/2,  =  0  (Dissociative,  Bates)  displays  behavior  inter¬ 

mediate  between  those  of  Fig.  4  and  6,  as  expected.  Also  shown  on  each  of  these 
maps  is  a  small  triangle,  the  hypotenuse  of  which  defines  contours  (straight  lines 


^An  improved  approximate  method  based  on  somewhat  different  arguments  is  also  dis¬ 
cussed  in  Appendix  B. 
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on  log- log  paper^)  of  constant  electron  mole  fraction.  As  will  be  shown*  a  know¬ 
ledge  of  this  straight  line  slope  allows  the  variation  of  electrical  conductivity 
in  the  frozen  section  to  be  followed  graphically  (see  Eq.  (17))  until  the  exit 
section  is  reached.  All  of  the  results  given  here  pertain  to  a  particular  family 
of  axisymmetric  nozzle  shapes  (frustum  of  cone  tangent  to  inside  of  torus  as  shown 
in  meridional  view  in  Fig.  7)  considered  to  be  a  fair  representation  of  those  of¬ 
ten  used  in  rocket  motor  practice.  The  divergence  half  angle  of  the  conical  sec¬ 
tion  was  assumed  to  be  15°  and  the  radius  of  curvature  at  the  throat  was  taken  to 
be  equal  to  the  throat  diameter. 

Figs.  4»  5  and  6  can  therefore  be  used  to  estimate  whether,  in  practice, 
the  recombination  parameter  C  is  small  enough  to  cause  significant  deviations 
from  ionization  equilibrium  within  the  nozzle  and,  if  so,  what  the  effect  on  the 
electrical  conductivity  at  the  exit  section  is  likely  to  be.  The  graphical  pro¬ 
cedure  is  illustrated  schematically  in  Fig.  8.  If  freezing  indeed  occurs  within 
the  nozzle  for  the  desired  pressure  ratio  and  estimated  recombination  parameter, 
then  the  appropriate  equilibrium  line  is  followed  to  the  point  where  it  intersects 

the  contour  of  constant  C  .  Thereafter,  one  follows  the  contour  of  constant 

r 

electron  mole  fraction  which  passes  through  this  point  until  the  desired  pressure 
ratio  is  reached.  The  enhancement  ratio,  cJ  /d  ,  due  to  ionization  nonequili- 

GX  GX  9  GC] 

brium,  is  then  calculated  from  the  ratio,  so  obtained  to  the  value 

*^ex,eq'^^c’  simply  following  the  equilibrium  line  down  to  the  desired 

pressure  ratio. 

3.  Effect  of  Physical  Size  (or  Thrust  Level)  on  the  Magnitude  of  the  Magneto- 
gasdynamic  Interaction  Parameter 

Quite  apart  from  the  scale  effect  introduced  by  ionization  nonequi librium 
phenomena,  the  effect  of  physical  scale  on  the  anticipated  value  of  the  magneto- 
gasdynamic  interaction  parameter  Q  at  the  nozzle  exit  is  explicitly  represented 
by  the  presence  of  the  ratio  in  Eq.  (9).  If  F^  is  the  thrust  level  of  the 


^When  the  electron  mole  fraction  is  constant,  Eqs.  (17)  and  (8)  combine  to  give 
d[ln  5]/d[ ln(l/p)]  =  (Y  ■  OA].  Thus,  contours  of  constant  electron  mole 
fraction  are  straight  lines  whose  slope  is  the  same  everywhere  on  the 
log  5  -  log  (l/p)  plane,  as  indicated  by  the  triangles  shown  in  Figs.  4,  5  and  6. 
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motor  in  question  (expressed  in  pounds)  and  p^ 
then  (in  accord  with  our  previous  choice  d 


for 


dt/d 


ref 


ref 


is 


is  the  chamber  pressure  in  psia 
-  1  meter)  a  convenient  relation 


> 


=  2.54  X  10-* 


*  i  F  i 


(20) 


i.e. t  for  a  fixed  propellant  combination,  chamber  pressure  and  expansion  ratio  we 
have  the  square  root  dependence:  ^.t/^ref 

It  should  be  remarked  that  while  values  of  Qj,gf  in  effect  compare  directly 
to  motors  of  the  same  geometric  size,  the  sequence  should  not  be  noticeably  al¬ 
tered  if  the  comparison  were  made  for  motors  of  equal  thrust  at  the  same  chamber 
pressure  level.  This  is  a  consequence  of  the  fact  that  (c*/c)^  exhibits  only  a 
very  sli^t  variation  in  going  from  one  propellant  combination  to  the  next.  For 
example,  (cVc)*^  only  varies  by  1.2  percent  for  the  first  six  propellant  combina¬ 
tions  listed  in  Table  III. 

4.  Effect  of  Nozzle  Expansion  Conditions  on  the  Magnitude  of  the  Magnetogas- 
dynamic  Interaction  Parameter  at  the  Nozzle  Exit  Section 

It  has  already  been  shown  that  the  value  of  the  interaction  parameter  Q  at 
the  nozzle  exit  section  can  be  computed  from 


Thus,  the  function  E  may  be  considered  to  account  for  all  of  the  remaining  effects 
of  the  gasdynamic  conditions  of  expansion  on  the  level  of  the  interaction  para¬ 
meter  Q.  This  composite  expansion  factor  has  been  plotted  once  and  for  all  for 
two  values  of  the  specific  heat  ratio  Y  =  ^p/^v  ~  1.20  and  1.30)  and  several 

values  of  the  ionization  parameter  ©A^  covering  the  pressure  ratio  range  of 
practical  interest.  Results  are  presented  in  Figs,  9  and  10. 
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V.  NUMERICAL  EXAMPLE;  TRANSVERSE  THRUST  FOR  A  300  PSIA  CHAMBER  PRESSURE, 
10,000  POUND  THRUST  OXYHYDROGEN  MOTORt  0.1^  POTASSIUM  SEED 


The  use  of  the  graphs  and  tables  already  presented  may  best  be  described 
with  the  help  of  a  numerical  example.  With  this  in  mind,  we  consider  the  combus¬ 
tion  chamber  specified  above  (in  the  section  title)  and  take  the  ambient  pressure 
to  be  one  atmosphere. 

1.  Deteimi nation  of 


For  a  one  percent  cesium  seed  is  obtained  from  Table  III  as  4.7  x  10“* 

For  a  tenth  of  a  percent  cesium  seed  the  result  would  be  4.7  x  10'  ■»A/io  or 
1.5  X  10“*  (see  Eq.  (4)).  To  convert  this  result  to  potassium,  we  make  use  of  Eq. 
(14)  noting  that  0^  -  0^^  =  50,110  -  44,900  =  5210°K  (Table  I),  T^  =  3437°K 
(Table  III)  and  =  132.91/39.10  =  3.4  (Table  I).  Therefore 

Or.f,K  =aref.C.  '  i(5210)/(3437)] 

or  o  1.3  r  I0r> 

2.  Effect  of  Electron-Ion  Recombination  Lag  (This  discussion  is  best  postponed 
until  after  Part  5) 

3.  Effect  of  Magnetic  Field  Strength 

Let  us  assume  that  a  magnetic  field  strength  of  3  Weber/meter*  is  available 
Then  (B/B^g^)*  =  (3/0*  =9  (see  Eq.(9)). 

4.  Effect  of  Physical  Scale  (or  Thrust  Level) 


From  Table  III  we  find  c*  =  2077  meter  (sec)"’  and  c  =  2920  meter  (sec)"’. 

Therefore,  since  the  thrust  level  F  is  10,000  lbs.  and  p  =  300  psia,  Eq.  (20) 

o  c 

gives 


ref 


=  2.54  X  10"*  (|)  (|^) 


,10,000,^  ^ 

^“"300"^ 
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5.  Determination  of  the  Equilibrium  Expansion  Factor»  E 

The  determination  of  E  involves  a  knowledge  of  0/T  ,  the  specific  heat 

C  Q 

ratio,  Y*  3nd  the  pressure  ratio  p  /p.  For  potassium  0  =  50,110  K  (Table  I), 

c 

so  that  0/T  =  50,110/3437  =  14.6.  The  pressure  ratio  is  given  by  300/14.7  = 

c 

20.4  (i.e. .  the  chamber  pressure  expressed  in  atmospheres).  Table  III  also  gives 
Y  =  Cp/Cy  so  that  Fig.  9  is  appropriate.  Introducing  these  values,  one 

finds  E  =  0.44.  With  the  exception  of  the  nonequilibrium  effect,  these  values  of 

snd  E  allow  an  estimate  to  be  made  of  the  interaction 
parameter  Q  to  be  expected  at  the  nozzle  exit  section.  Since  ^  1 

we  can  then  immediately  state  the  result 

or  Q  (1.3  r  10-*)(9)(0.140)(O.U)  =0.72  x  lO'* 

The  extent  to  which  nonequilibrium  ionization  will  increase  the  actual  value  of 
the  interaction  parameter  Q  above  its  minimum  value  of  0,72  x  10"*  will  now  be 
considered  in  the  light  of  the  analysis  of  Appendix  B, 

2,  (Continued)  Effect  of  Electron-Ion  Recombination  Lag 


To  assess  the  importance  of  electron-ion  recombination  lag  it  is  necessary 
to  estimate  the  most  important  electron  loss  mechanism**  and,  on  the  basis  of  this 
mechanism,  evaluate  the  most  probable  value  of  the  recombination  parameter  de¬ 
fined  by  Eq,  (l8).  Before  doing  this,  however,  it  can  be  verified  that  for  a  giv¬ 
en  chamber  pressure,  seed  material  and  concentration,  small  values  of  will  be 
favored  by  propellant  combinations  with  low  values  of  the  characteristic  velocity 
c*.  This  is  because  the  effect  of  reduced  chamber  temperature  on  the  electron 

number  density  n  in  the  chamber  outweighs  the  effect  of  temperature  on  the  re- 
e  f  c 

combination  coefficient  a  or  the  direct  effect  of  reduced  c*  on  C_.  It  is  also 

r  ,c  r 

clear  that  small  values  of  the  recombination  parameter  will  be  favored  by  small 
(throat  diameter)  motors.  What  is  least  clear,  unfortunately,  is  the  establish¬ 
ment  of  the  most  probable  value  of  the  recombination  coefficient  (a,  _)  in  the 

TfC 

chamber,  and  the  law  governing  its  temperature  and  pressure  variation  (i .e, ,  the 

coefficients  w  and  u  ). 

P  T 
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The  solid  curve  shown  in  Fig.  11  gives  the  pressure  dependence  of  the  ion 
recombination  coefficient  for  air  at  ambient  temperature.^  The  curve  may  be  seen 
to  have  two  distinct  branches  separated  by  a  transition  region.  The  low  pressure 
branchy  labeled  Thompsony  300°K>  is  believed  to  scale  according  to  ywhereas 

the  high  pressure  branch  labeled  Langevin,  300°Ky  is  believed  to  scale  according 
to  T’/*.  The  displaced  lines  marked  3000°K  then  represent  the  corresponding 
branch  positions  at  a  temperature  level  more  representative  of  rocket  motor  prac¬ 
tice.  If  the  a^(p;T)  curves  for  other  gases  retain  this  same  shape  at  elevated 
temperaturesy  then  the  Thompson  coefficient  should  prevail  at  rocket  motor  pres¬ 
sures  and  this  value  would  be  considerably  smaller  than  the  value  obtained  by 
extrapolating  the  Langevin  branch.  Moreovery  since  pressure  and  temperature  will 
be  related  as  in  Eq.  (8)y  the  Thompson  branch  should  prevail  throughout  the  noz¬ 
zle.  This  is  illustrated  by  the  locus  of  recombination  coefficients  a  (p)y  super- 

r 

imposed  on  Fig.  11,  and  calculated  for  isentropic  expansion  from  a  300  psia, 

3000®K  chamber  with  Y  ~  The  implication  is  thaty  of  the  sets  of  coefficients 

w^yWp  introduced  earliery  the  set  which  gives  rise  to  the  earliest  departures  from 
equilibrium  is  indeed  applicable  throughout  the  nozzle.  It  can  therefore  be  said 
that  predictions  based  on  the  Thompson  mechanismy  by  giving  the  smallest  absolute 
estimates  of  the  recombination  rate  parameter  and  the  largest  anticipated  de¬ 
partures  from  equilibrium  (compare  Figs.  4y5y6)y  represent  the  most  optimistic 
case  for  magnetogasdynamic  applications,  i.e. ,  lead  to  the  largest  expected  values 
of  the  nonequilibrium  enhancement  .  Returning  to  the  particular  numer- 

©X  ©X  j  ©C| 

ical  example  treated  in  this  section,  it  can  now  be  shown  that  even  in  the  most 
optimistic  case  no  significant  nonequilibrium  enhancement  of  the  electrical  con¬ 
ductivity  at  the  exit  section  should  be  expected.  This  can  be  demonstrated  as 
follows. 

The  reference  electron  number  density  n  given  in  Table  III  is  0.74  x  lO’* 

©  j  c 

electrons/cm*.  When  corrected  for  0. 1/f  potassium  seed,  this  becomes  about  2  x  lO’^ 
electrons/cm*.  The  recombination  coefficient,  a  ,■ under  these  pressure-temper- 

Tj  C 

ature  conditions  is  expected  to  be  of  the  order  of  lO”*  cm^(sec)“’  (see  Fig.  ll). 

Therefore,  the  characteristic  time  for  electron  decay  will  be  = 

(a  n  -  5  X  10“*  sec.  On  the  other  hand,  the  characteristic  flow  time 

ryC  e,c 

Tfio^  =  d^/c*  is  found  to  be  0.140/2077  =  0.67  x  10"*  sec.  Therefore,  the  recom¬ 
bination  parameter,  =  ’^f low^"''chem’  order  of  1.3  x  10*.  Consulting 


^This  graph  is  extracted  from  Fig.  77,  p.  135  of  Von  Engel**,  and  strictly  applies 
to  air  at  room  temperature  (see  Appendix  B). 
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Fig.  6,  we  notice  that  such  a  contour  would  cross  the  (interpolated)  equilibrium 

curve  corresponding  to  0/T  =  14.6  at  a  value  of  p  /p  approximately  equal  to 

c  c 

18.5  (i .e. «  fairly  near  the  exit  section  value  of  20.4).  Therefore,  even  in  the 
most  optimistic  case  for  this  numerical  example,  departures  from  ionization  equi¬ 
librium  at  the  nozzle  exit  section  are  expected  to  be  small,  being  of  the  order 
of  only  16  percent.  We  conclude  that  to  within  the  accuracy  of  the  present  anal¬ 
ysis,  the  inequality  sign  appearing  in  the  earlier  statement:  Q  >y0.72  x  10"* 
could  be  removed.  At  best  we  may  take  Q  =  0.84  x  10"*. 

It  is  now  possible  to  estimate  the  magnetogasdynamic  control  force  that 
would  result  from  the  application  of  a  30,000  gauss  (3  Weber/meter*)  magnetic 
field  (inclined  at  45°)  to  the  exhaust  jet.  If  the  axial  extent  of  the  field  is 

comparable  to  the  jet  diameter,  then  L/d  ~  1  and  F, /F  is  found  from  Eq.  (A18)  of 

t  X  o 

Appendix  A,  to  be  approximately  Q  =  4.2  x  10“^.  Thus  0.42  percent  of  the 
gross  thrust  would  be  available  for  control  purposes.  Since  the  undisturbed 
thrust  level,  F  ,  is  10,000  pounds,  the  anticipated  control  force,  Fi  ,  would  then 

O 

be  approximately  42  pounds  for  this  case.  A  field  strength  of  lOOiCXXD  gauss  (lO 
Weber/meter^)  would  increase  this  value  to  (l0/3)*  •  (42)  or  about  470  pounds»  i 
slightly  less  than  5  percent  of  the  gross  thrust  of  the  engine. 

VI.  DISCUSSION 


An  examination  of  the  equations  and  figures  presented  in  previous  sections 
allows  some  general  conclusions  to  be  drawn  about  the  effects  of  various  para¬ 
meters  on  the  anticipated  magnetogasdynamic  control  force.  While  it  is  quite  ob¬ 
vious  that  large  magnetic  field  strengths  are  desirable,  the  effects  of  such  para¬ 
meters  as  physical  scale,  chamber  pressure,  specific  heat  ratio  etc.,  are  not 
quite  so  evident,  particularly  when  the  possibility  of  nonequilibrium  ionization 
at  the  exit  section  is  allowed  for.  In  accord  with  the  previous  formulation,  the 
individual  influence  of  these  parameters  (other  than  field  strength)  on  the  inter¬ 
action  parameter  Q  can  be  discussed  by  forming  the  ratio 


^  ^  V  i  ^V^ex,eq)’ 

Q  -Qref  ■  ^  ‘  (<^ex/^ex,eq) 


(21) 


^If  Q  had  exceeded  10"*  the  use  of  Fig.  2  would  be  more  convenient. 
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where  the  primed  variables  may  be  taken  to  represent  values  of  the  individual 
quantities  after  a  particular  parameter  (say,  0/T  ,  Y  *  Pc^^ex’  '^t'^'^ref^  been 
changed.  Since  the  control  forces  are  monotonic  in  the  interaction  parameter  Q, 
any  factor  which  (all  things  considered)  causes  Q'  to  exceed  Q  would  be  expected, 
in  practice,  to  favor  larger  MHD  control  forces. 

1.  Effect  of  Q/T  on  Interaction  Parameter 

c 

It  is  appropriate  that  the  effect  of  this  parameter  be  discussed  first 
since  it  is  found  to  dominate  all  others.  Indeed,  a  consideration  of  its  magni¬ 
tude  accounts  for  the  effect  of  the  remaining  parameters  on  the  interaction  para¬ 
meter  Q. 

If  the  only  difference  between  two  rocket  motors  was  the  value  of  the 

characteristic  ionization  temperature  0  corresponding  to  the  seed  material,  then 

the  values  of  Qj,gf»  E  and  ^ex/'^ex,eq  likewise  be  different.  and  E 

would  be  larger  for  the  low©  seed  material  motor  since  Eq.  (12)  shows  that 

Q  ,  oc  exp  (-  "i- 0/T  ),  while  Figs.  9  and  10  reflect  the  strong  effect  of  0/T 
r©T  c  c 

on  the  factor  (^gx^^c^eq  appearing  in  the  equation  for  E.  However,  it  cannot  be 
concluded  that  Q  would  be  larger  until  the  nonequilibrium  factor  «/,) 

6X  cX  y  ©q 

investigated.  This  is  particularly  true  since  low  values  of  0  tend  to  increase 
the  magnitude  of  the  recombination  parameter  C^,  thereby  suggesting  smaller  de¬ 
partures  from  equilibrium.  In  addition,  the  value  of  the  normalized  equilibrium 
electrical  conductivity  (cJ  /d  )  will  be  larger,  which  also  suggests  smaller 

@X  C 

values  of  the  nonequilibrium  enhancement  C  ^  When  these  factors  are  com- 

©X  ©Xf©C) 

binedy  hovyevety  it  is  found  that  the  equilibrium  benefits  of  low  characteristic 
temperature  ©  more  than  compensate  for  the  nonequilibrium  gains  associated  with 
larger  values  of  0  for  the  parameter  range  likely  to  be  of  interest  in  chemical 
propulsion.  Thusy  even  in  the  nonequilibrium  regime  it  would  appear  to  be  advan¬ 
tageous  to  use  the  seed  material  with  the  lowest  possible  ionization  potential. 

2.  Effect  of  Specific  Heat  Ratio 

For  motors  in  which  equilibrium  ionization  is  likely  to  be  maintained 
throughout  the  nozzley  the  strong  effect  of  specific  heat  ratio  which  shows  up  in 
Figs.  9  and  10  can  be  traced  to  the  steep  temperature  dependence  of  the  electri¬ 
cal  conductivity  for  realistic  values  of  the  parameter  0/T^.  To  illustrate  this 
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point,  let  ut  consider  the  effective  value  of  the  exponent  n  in  the  power  law^"^ 

CJ  oc  t"  (22) 

assuming  that  cesium  (i.e. .  the  element  having  the  lowest  available  value  of  0) 
is  used  as  the  seed  material.  For  the  purpose  of  comparing  effective  values  of 
n  likely  to  be  realized  in  rocket  motor  practice,  it  is  convenient  to  consider 
a  temperature  (drop)  ratio  which  is  characteristic  of  each  propellant  gas  combin¬ 
ation,  say  2/(Y  +  l),  which  will  be  recognized  at  the  ratio  of  the  throat  temper¬ 
ature  to  the  chamber  temperature  for  an  Isentropic  expansion.  This  temperature 
drop,  if  unaccompanied  by  a  partially  compensating  pressure  drop,  would  give  rise 
to  values  of  n  which  can  be  calculated  from  a  relation  easily  derived  from  Eq.  (7) 


JL  .  Hyh-1)  -  1  2 

^^c  ln[i(Y+l)]  4 


(23) 


Values  of  the  index  n  were  calculated  for  the  first  six  propellant  combinations 
listed  in  Table  III  and  with  the  results  given  in  Table  VI. 


TABLE  VI 

EFFECTIVE  VALUES  OF  d(ln  0)/d(ln  T)  FOR  CESIUM 
SEEDING  IN  ROCKET  MOTOR  APPLICATIONS 


Propellant 

T 

c 

Y 

n 

L0X-NH3 

3010 

14.8 

1.23 

8.6 

Lox-Gasolene 

3295 

13.6 

1.22 

7.9 

Lox -Hydrazine 

2750 

16.3 

1.23 

9.4 

Lox-Kerosene 

3350 

13.3 

1.24 

7.8 

Fj-Hydrazlne 

4430 

10.1 

1.33 

6.2 

N^O^-Hydrazine 

2975 

15.0 

1.27 

8.8 

Thus  the  equilibrium  electrical  conductivity  <J  would  vary  as  almost  the  10th 
power  of  the  absolute  temperature  for  lox-hydrazine  and  as  about  the  6th  power  of 
the  temperature  for  f lourine-hydrazine.  For  seed  materials  other  than  cesium. 
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these  exponents  would  be  even  larger  (almost  in  the  ratio 

In  relating  this  behavior  to  the  anticipated  effect  of  the  specific  heat 
ratio  Y  on  the  interaction  parameter  Q,  it  is  then  not  surprising  that  values  of 
Y  giving  the  smallest  temperature  change  for  a  given  pressure  drop  are  advanta- 
geous.  Since  Eq.  (8)  indicates  that  T  oc  p''  ”,  it  can  be  seen  that  small 

values  of  the  specific  heat  ratio  Y  give  the  largest  values  of  the  equilibrium 
expansion  factor  E.  Hence,  all  other  factors  being  equal,  the  interaction  para¬ 
meter  Q  at  the  exit  section  will  be  maximized  by  propellant  combinations  having 
small  effective  heat  capacity  ratios. 

3.  Effect  of  Expansion  Pressure  Ratio  for  a  Class  of  Equal  Thrust  Motors 
Expanding  to  the  Same  Back  Pressure 

Increased  chamber  pressure  will  decrease  the  factors  d,  and  E  (in  the 

0A_  range  of  practical  interest)  and,  as  such,  tend  to  reduce  the  strength  of  the 
c 

magnetogasdynamic  interaction  at  the  exit  section.  This  is  partially  offset  by 
the  fact  that  higher  pressures  give  rise  to  larger  values  of  the  nonequilibrium 
enhancement  ,  despite  the  fact  that  (according  to  the  Thompson  mech- 

anism)  the  recombination  paramoter  will  also  incraase.  Again,  tha  decrease 
in  the  equilibrium  performance  parameters  will  outweigh  the  nonequilibrium  gain 
with  the  result  that  reduced  chamber  pressures  will  generally  favor  increased  mag¬ 
netic  interaction  for  motors  of  the  same  thrust  level.  This  conclusion  is  also 
valid  in  comparing  motors  of  the  same  physical  size  (for  which  a  chamber  pressure 
increase  implies  a  thrust  increase), 

4.  Scale  Effect  for  Motors  of  the  Same  Pressure  Level 


The  vertical  spacing  of  the  nonequilibrium  contours  (lines  of  constant  C^) 
shown  in  Fig.  6  suggests  that,  even  in  the  nonequilibrium  regime,  increased  size 
will  always  favor  increased  magnetic  interaction.  This  is  a  consequence  of  the 
fact  that  reducing  the  throat  diameter  and,  hence  by  an  order  of  magnitude, 
does  not  cause  enough  of  an  increase  in  the  exit  value  of  the  conductivity  to  out¬ 
weigh  the  explicit  cJ/d^  dependence  appearing  in  Eq.  (9).  In  the  equilibrium  re¬ 
gime,  this  latter  dependence  on  physical  size  is  indeed  the  only  one. 

Summarizing  the  results  of  each  of  paragraphs  (l)  through  (4)  above,  it  may 
be  said  that,  in  the  range  of  rocket  motor  interest,  the  following  characteristics 
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appear  to  be  generally  desirable  for  MHD  thrust  vectoring  applications: 

i  small  values  of  ©A^  (i.e. ,  low  ionization  potential  seed 
materials  in  "hot"  chambers) 
ii  small  values  of  the  specific  heat  ratio  Y 
iii  moderate  rather  than  high  chamber  pressures 
iv  large  physical  size 

5.  Accuracy  of  Physical  Model 

With  regard  to  the  absolute  values  predicted  by  the  use  of  the  analysis* 
several  comments  should  be  made  about  the  accuracy  of  the  physical  model.  In  gen¬ 
eral,  most  departures  from  the  model  that  might  occur  in  practice  would  tend  to 
reduce  the  magnitude  of  the  obtainable  control  force.  Factors  of  this  type  include: 

(a)  magnetic  field  inhomogeneity 

(b)  ion  slip  and  Hall  effects** 

(c)  electron  attachment**  to  form  negative  ions 

(d)  reactions  of  seed  material  with  combustion  products 

(e)  non-uniformity  of  profiles  of  mass  velocity  and  electrical  conductiv¬ 
ity  at  the  nozzle  exit 

(f)  electron  removal  from  the  jet  by  diffusion 

Apart  from  factors  of  this  type,  the  accuracy  of  the  predictions  may  be  limited 
for  particular  systems  for  such  reasons  as: 

(g)  the  expansion  process  cannot  be  well  represented  by  any  single  value 
of  the  isentropic  exponent  Y 

(h)  the  temperature  dependence  of  the  electron-neutral  cross-section 
is  not  negligible*' 

(i)  the  pressure  and  temperature  dependence  of  the  recombination  coeffi¬ 
cient  for  the  effective  electron  removal  process  are  not  well  des¬ 
cribed  by  a  single  power  law 

(j)  the  absolute  value  of  the  recombination  coefficient  under  chamber  con¬ 
ditions  is  not  well  known 

etc. 


Any  power  law  dependence  would  be  easily  incorporated  into  the  present  analysis. 
For  the  electron  interaction  with  water  vapor*,  Q  ~  T"’  would  be  more  accurate 
than  the  law  used  in  constructing  Figs.  4,5,6. 
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Finally,  it  should  be  remembered  that  even  if  the  kinetics  of  electron  removal 
were  accurately  known,  most  of  the  nonequilibrium  results  presented  here  (Fig.  4» 
5»6)  are  approximate.  Where  high  accuracy  is  essential,  exact  solutions  to  Eq. 
(B-8)  should  be  used.  It  is  felt  that  the  results  given  here  (see  also  Appendix 
B)  correctly  parameterize  the  problem,  correctly  represent  the  dominant  trends, 
and  enable  rapid  estimates  to  be  made.  It  should  also  be  borne  in  mind  that  even 
the  predictions  of  exact  solutions  (see  Appendix  B)  should  not  be  applicable  for 
nozzle  shapes  having  area  change  histories  markedly  different  from  that  considered 
here  (see  Fig.  lO). 


VII.  CONCLUDING  REMARKS 


Recapitulating,  several  universal  curves  have  been  presented  in  Sections 
I,  II  and  III  to  facilitate  the  prediction  of  control  forces  that  can  be  achieved 
by  magnetohydrodynamic  means  in  rocket  motor  practice.  A  set  of  equilibrium  cal¬ 
culations  also  required  for  this  purpose  has  been  presented  for  21  propellant  com¬ 
binations  assumed  to  be  "seeded"  with  small  amounts  of  an  easily  ionized  material. 
Since,  in  principle,  the  nozzle  expansion  process  can  result  in  electrical  conduc¬ 
tivities  at  the  exit  section  considerably  greater  than  those  corresponding  to  local 
thermodynamic  equilibrium,  a  graphical  method  has  been  introduced  to  enable  quan¬ 
titative  estimates  to  be  made  of  the  extent  of  this  electron  recombination  lag  as 
well  as  the  conditions  under  which  it  can  (very  likely)  be  ignored.  Combining 
these  results  as  illustrated  in  Sections  V  and  VI,  it  has  been  possible  to  extract 
some  general  information  as  to  the  independent  factors  favoring  maximum  exhaust 
jet  interaction  with  the  applied  magnetic  field.  Some  of  the  major  factors 
(other  than  large  field  strength)  giving  rise  to  favorable  conditions  are: 

(a)  the  presence  of  low  ionization  potential  materials  in  "hot"  chambers  (low 
0/T^),  (b)  small  values  of  the  effective  specific  heat  ratio  Yt  (c)  moderate  rather 
than  extremely  large  chamber  pressures,  and  (d)  large  physical  size.  It  appears 
from  the  calculations  carried  out  thus  far  that  control  forces  of  the  order  of  0.1 
to  1.0  percent  of  the  gross  thrust  are  attainable,  by  applying  large  magnetic 
fields  (l  to  3  Weber/meter*)  to  the  exhaust  jets  of  rocket  motors  in  the  10,000 
pound  thrust  class  over  a  length  comparable  to  the  exit  diameter.  Upper  limits 
to  the  seed  concentrations  are  considered  to  be  in  the  range  of  0.1  to  1.0  per¬ 
cent.  With  the  use  of  superconducting  magnet*  coils  (perhaps  capable  of  generat¬ 
ing  100,000  gauss  fields),  side  thrusts  of  the  order  of  several  percent  of  the 
gross  thrust  are  predicted.  Interestingly  enough,  nonequilibrium  effects  are 
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tentatively  estimated  to  be  quite  small  over  a  wide  range  of  conditions  likely  to 
be  encountered  in  rocket  motor  practice  and,  moreover,  the  factors  favoring  large 
nonequilibrium  gains  are  predominantly  those  associated  with  small  equilibrium 
values  of  the  interaction  parameter.  The  net  effect  is  that  control  force  expec¬ 
tations  of  this  order  of  magnitude  generally  apply  to  both  the  equilibrium  as  well 
as  nonequilibrium  cases  envisioned.  Whether  these  control  forces  are  sufficient 
to  justify  a  serious  effort  in  this  direction  remains  to  be  seen,  since  economic 
and  technological  factors  are  involved  which  have  not  even  been  touched  upon  here. 
Insofar  as  the  material  presented  here  would  provide  perspective  and  input  infor¬ 
mation  for  such  a  study,  it  is  hoped  that  this  report  will  prove  useful. 
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APPENDIX  A 

TRAJECTORY  ANALYSIS;  MAGNETOGASDYNAMIC  JET  DEFLECTION 

Using  an  approach  analogous  to  that  adopted  by  Dimmock*’  etal.  for  finding 
the  deflection  of  flame  gas  in  the  presence  of  crossed  electrical  and  magnetic 
fields,  it  is  possible  to  obtain  estimates  for  the  extent  of  jet  deflection  in 
the  presence  of  ja  magnetic  field  alone.  In  this  case,  we  assume  that  a  fluid  par¬ 
cel  of  unit  volume  and  mass  p  will  be  acted  upon  solely  by  the  Lorenz  force  j  x  §, 
where  the  current  3  is  induced  by  the  gas  motion  through  the  field  in  accord  with 
3  =  CJ(c  X  B).  If  the  magnetic  field  is  in  the  x-y  plane,  then  the  equations  of 
motion  for  the  fluid  parcel  shown  in  Fig.  1  take  the  form 

pV  =  -  ct(uBy  -  vB^)By  (A-1) 

py  =  cJ(uBy  -  vB^)B^  (A-2) 

The  initial  conditions,  corresponding  to  the  undisturbed  gas  flow^are  taken  to  be 

x(0)  =  0  x(o)  =  c  (A-3) 

o 

y(0)  =  0  y(0)  =  0  (A-4) 

Before  proceeding,  it  is  convenient  to  introduce  the  following  definitions  and 
nondimensionali ration 


X  E  x/d  TEC  t/d 

°  (A-5) 

Y  =  y/d  Q  =  cJB*d/(pc^) 

We  further  take  the  impressed  magnetic  field  to  have  a  magnitude  B  and  an  inclin¬ 
ation  angle  0  with  the  x  axis.  Then,  using  primes  to  represent  differentiation 
with  respect  to  x,  the  governing  set  of  equations  assumes  the  form 

Q-’  X"  =  -  X'  sin*0  +  Y'  sin  9  cos  0  (A-6) 

Q”’  Y"  =  X'  sin  0  cos  9  -  Y'  cos*9  (A-7) 
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with  the  corresponding  initial  conditions 


X(0)  =0  X’  (0)  =  1 

(A-8) 

Y(0)  =0  Y'  (0)  =  0 

(A-9) 

In  writing  Eqs,  (a-6)  and  (A-7),  it  has  been  assumed  that  the  induced  magnetic 
field  may  be  neglected.  The  solution  to  Eqs.  (A-6)  and  (A-7)»  subject  to  condi¬ 
tions  (A-8)  and  (A-9),  is  found  to  be 


X(t)  =  Q"’  sin*0[l  -  exp(-  Qr)]  +  T  cos*0  (A-IO) 

X’(t)  =  1  -  sin*0[l  -  exp(-  Qr)]  (A-ll) 

y(t)  =Q"'’  sin  e  cos  6[Qt  -  [l  -  exp(-  Qt)]J  (A-12) 

Y’(t)  =  sin  0  cos  0[l  -  exp(-  Qr)]  (A-13) 


Since  the  jet  turning  angle  p  will  be  given  by  tan"’[Y' (t)/X' (t)]  and  the  result¬ 
ing  thrusts  are  given  respectively  by  =  Y*  (-t)  and  =  X’(x)»  we  have 


p  =  tan" 


,rsin  0  cos  0[l  -  exp(-  Qx)]"^ 


^rsin  B  cos  HLl 
^  1  -  sin*0[l  - 


0[l  -  exp(-  Qt)] 


(A-14) 


Fy/F^  =  sin  0  cos  9[l  -  exp(-  Qr)]  (A-i5) 

F^/F^  =  1  -  sin*e[l  -  exp(-  Qr)]  (A- 16) 

Finally,  from  Eq.  (A-ll)  the  time  spent  in  a  field  of  axial  extent  L  is  found  to 
be  the  solution  of 


L/d  =  Q’’  sin*6[l  -  exp(-  Qx)]  +  x  cos*0  (A-17) 

For  small  Qx  f i . e .  sma  1 1  0  B*x/(pc^)]  x  and  L/d  become  indistinguishable  from  one 

another  and  F  /F  (written  hereafter  as  F  /F  )  approaches 
y  o  X.  0 
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Fj^/F^  ->  sin  6  cos  0  •  (L/d)  •  Q  (A-18) 

The  relation  between  the  transverse  thrust  F,  /F  ,  interaction  parameter  Q,  and 

field  extent  parameter  L/d  is  shown  graphically  in  Fig.  2  for  the  optimum  field 

inclination*  6  =  45°.  For  values  of  Q  smaller  than  10"**  the  asymptotic  result 

(A-18)  can  be  used  with  good  accuracy,  and  this  is  usually  the  case  in  practice. 

Note  that  in  this  extreme  the  transverse  thrust  will  increase  linearly  with  Q, 

and  that  if  0  =  45°  and  L/d  =  1,  then  Fj/F^  i  Q.  The  asymptotic  extreme 

Q -*  oo  is  of  more  theoretical  interest.  In  this  extreme  Fi/i^  -*^,  corresponding 

J.  0 

to  a  maximum  jet  deflection  equal  to  the  45  field  inclination. 
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APPENDIX  B 

ELECTRON-ION  RECOMBINATION  LAGS 


Formulation 

To  predict  the  extent  of  departures  from  equilibrium  in  nozzle  flows* 
one  must  examine  the  conservation  equation  governing  changes  in  the  ion  or  elec¬ 
tron  concentration.  Since  we  are  really  interested  in  the  electrical  conducti¬ 
vity  of  the  gas  and  its  variation  through  the  nozzle,  it  is  convenient  to  trans¬ 
form  this  equation  into  one  in  which  the  conductivity  cJ  is  the  dependent  variable, 
in  accord  with  Eq.  (2).  The  particular  case  to  be  considered  here  is  one  in  which 
a  seed  material,  present  in  dilute  amounts,  contributes  most  of  the  electrons  to 
the  gas  via  the  reaction 

S  -*  s'*"  +  e  (B-l) 

In  the  course  of  the  expansion,  the  reverse  reaction  begins  to  dominate  (B-l) 

s’*"  +  e  ->■  S  (B-2) 

and  if  the  expansion  is  rapid  enough,  S^,  e,  and  S  may  not  be  in  equilibrium  with 
one  another  downstream  of  the  coinbustion  chamber.  Regardless  of  the  actual  molec- 
ularity,  the  kinetics  of  reaction  (B-2)  are  conventionally  described  by  a  phen¬ 
omenological  law  of  the  form 


dn 
_ e 

=  -  a  n  n. 
dt  r  e  + 


(B-3) 


where  is  referred  to  as  the  recombination  coefficient.  With  the  condition  of 
electrical  neutrality  (n^  =  n^),  Eq.  (B-3)  becomes 


dn 
_ e 

dt 


(B-4) 


It  follows  that  at  any  point  in  the  nozzle  the  net  rate  of  electron  pro¬ 
duction  can  be  expressed  as  -  a  [ (n  )*  -  (n  _„)*]  in  units  of  number  of  elec- 

r  6  6  f  ©H  , 

irons  per  unit  time  and  volume.  Therefore,  the  mass  source  term  w  appearing  in 
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the  conservation  equation  for  electrons  is  of  the  form 


;  =-  o  m  r(n  -  (n^ 

e  r  e  e 


(B-5) 


where  m  is  the  mass  of  the  electron.  In  terms  of  the  mass  fraction,  m^  n  /p,  of 
0  Be 

electrons,  the  steady  flow  conservation  equation  therefore  becomes 


(B-6) 


As  already  suggested,  it  is  convenient  to  nondimensionalize  and  transform  Eq.  (B-6) 
into  one  governing  changes  in  the  electrical  conductivity  (3  itself.  To  do  this, 
we  introduce  the  following  definitions 


p  =  p/Pj 

f  =  tA^ 
5  =  cj/p^ 

X  =  a  /o 


r,c 


u  =  u/c* 

X  =  x/D^ 

^eq  “  ‘^eq^*^c,eq 
=  D  T 


(B-7) 


and  denote  the  group  a  n  D  /c*,  which  naturally  makes  its  appearance,  as 

r>  c  G  >  c  t 

the  recombination  parameter,  C^.  We  note  that  is  nothing  more  than  the  ratio 


of  the  characteristic  flow  time 


T-,  =  dJc* 

flow  r 


to  the  electron  half  life 


T  .  =  (q  n  in  the  chamber.  The  resulting  equation  for  the  nonequili- 

ch0in  f  9  ^  ^9^ 


brium  conductivity  distribution  3(x)  becomes 


§  =  -  c  (5)“"*’  (C)-’[(a)»  -  (3„)']  -  i  I  4  (B-8) 

dx  Cm  T  dx 

This  equation  is  in  a  convenient  form  since  the  variation  of  5^^  through  a  nozzle^ 
has  already  been  calculated  and  discussed  (see  Section  III).  Moreover,  the  rela¬ 
tions  between  u,  p,  and  T  are  particularly  simple  for  isentropic  nozzle  flow,  i.e. 

Y 

p  =  (B-9) 


and 


G  =  [2Y/(Y  -  (1  -  t)^ 


(B-10) 
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Y+t 

where  Q(Y)  i  (Y)^[2/(Y  +  (B-ll) 

is  a  well  tabulated  function^^  of  the  specific  heat  ratio  Y.  Exact  solutions  to 
Eq.  (B-8)  could  therefore  be  obtained  by  numerical  integration  for  any  class  of 
nozzle  shapes.  However,  before  embarking  on  such  a  program,  it  was  decided  to  ex¬ 
tract  some  approximate  results  without  performing  any  integrations.  Indeed,  the 
results  described  below  (which  are  reported  in  the  text)  were  obtained  using 
approximations  of  the  type  which  have  been  successfully  introduced  by  Bray  and 
Hall  etal.  for  chemically  reacting  nozzle  flows.  A  comparison  of  these  results 
with  a  more  exact  treatment  will  also  be  included, 

"Sudden  Freeze”  Approximation 

The  principle  underlying  these  methods  is  that  the  transition  region  be¬ 
tween  equilibrium  and  frozen  flow  is  frequently  small  enough  to  disregard.  The 
actual  distribution  of  species  in  the  nozzle  is  therefore  often  well  represented 
by  patching  the  equilibrium  and  frozen  solutions  together  at  a  suitable  "freeze 
point".  Methods  of  determining  this  freeze  point  differ  somewhat}  however,  for 
the  present  purposes,  any  of  the  methods  commonly  used  should  suffice.  What  is 
important  is  that  this  point,  and  the  corresponding  electron  mole  fraction,  can 
be  obtained  from  the  properties  of  the  equilibrium  curve  itself. 

A  method  of  determining  the  freeze  point  paralleling  that  of  Ref.  17  would 
proceed  as  follows.  At  the  freeze  point  it  is  approximately  true  that 


-  4  ~  C  (p)  P  (T)  (8  )2  +  i  4  (B-12) 

dx  ^  eQ  T  dx 


If,  in  that  region,  we  introduce  the  "average"  slope 


+  (?)  ] 

dx  dx  ^  dx  ^ 


(B-13) 


Then  freezing  should  occur  where  the  following  equality  is  achieved 
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iu(p)  (t)  (5 


eq 


)-[ 


d 

dx 


(B-U) 


For  a  nozzle  of  any  prescribed  shape,  the  quantity  on  the  left  hand  side  of  this 
equation  can  be  tabulated  as  a  function  of  Mach  number  or  p.  Thus  the  value  of 
p,  and  hence  the  corresponding  conductivity  5^,  can  be  obtai»ied  for  any  desired 
value  of  the  rate  parameter  C^, 

Calculation  Procedure 


To  carry  this  procedure  through  in  a  realistic  case,  calculations  were 
perfoonned  for  the  nozzle  shape  shown  in  Fig.  7.  This  axisymmetric  nozzle  includes 
a  conical  expansion  section  with  a  half  angle  of  15^.  The  cone  is  in  turn  tan¬ 
gent  to  the  inside  of  a  torus  which  comprises  the  entire  converging  section  (and 
a  small  portion  of  the  expansion  section).  The  torus  is  sized  such  that  the 
radius  of  curvature  at  the  throat  of  the  resulting  nozzle  is  equal  to  its  throat 

diameter  D  .  It  is  then  a  simple  matter  to  express  the  area  and  its  spatial  rate 
^  -  ,7r 

of  change  within  this  nozzle.  For  x  1  -  cos(^  +  the  dimensionless  area 
distribution  is  found  to  be 


A  H  A/A^  =  [1  +  2  (1  -  sin*>^]^ 


(B-15) 


where'll  and  x  are  related  in  accord  with  x  =  1  -  cos  ^  .  Therefore,  in  this 
range 


dx 


4 

tan'Z^ 


On  the  other  hand,  in  the  conical  section 


4  =  4  tan  <t  (A)^ 
dx 


(B-16) 


(B-17) 


The  derivatives  required  to  proceed  are  dO^^/dx  and  df/dx.  It  is  conven¬ 
ient  to  calculate  these  quantities  by  decomposing  them  in  accord  with  the  chain 
rule  as  follows 
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r|  1  df 

dx  ~Uaf  /-  dM  lap  /.  dMJjdMl 


(B-18) 


dx 


(^ll 

dM/l 

idM/ 

dx 


(B-19) 


An  expression  for  each  individual  derivative  is  now  available.  In  particular,  it 
can  be  verified  from  Eq.  (7)  that 


da 


eq 
df  j 


[I- iir  ^  I  Iff] 


8CJ. 


_ eq 

l9p  / 


=  ..1.5 

2p  «•=! 


(B-20) 


(B-21) 


From  the  isentropic  flow  relations,  one  can  obtain  the  following  three  derivatives 
with  respect  to  Mach  number 


g  =  -  rr  -  1)  (f)^  « 


(B-22) 


dM 


^  -  Y  p  T  M 


(B-23) 


I  =  (M=  -  1)  f  (i)  (B-24) 

Using  the  relations  given  above  and  a  set  of  isentropic  flow  tables^t  it 
is  then  possible  to  make  a  plot  of  the  left  hand  side  of  Eq.  (B-14)  as  a  function 
of  the  Mach  number  for  any  choice  of  the  coefficients  and  (i  .e.  <  for  any  re¬ 
combination  mechanism).  Having  done  this  (see  Fig.  12) »  the  conductivity 
corresponding  to  any  prescribed  value  of  the  rate  parameter  can  be  read  off  a 
figure  such  as  Fig.  6.  If  the  carrier  gas  expansion  proceeds  past  this  pointt 
then  d  will  thereafter  vary  in  accord  with  Eq.  (17)  given  in  the  text.  Because 
pressure  and  temperature  bear  a  power  law  relationship  with  one  anotheri  a 
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convenient  way  of  presenting  these  results  is  to  superpose  contours  of  constant 
recombination  parameter  on  a  log-log  plot  of  equilibrium  conductivity  versus 
expansion  pressure  ratio  (c.f.  Figs.  4»5>6).  If  the  contour  corresponding  to  the 
best  estimate  of  (for  a  particular  problem)  crosses  the  equilibrium  curve  be¬ 
fore  the  exit  pressure  is  reached,  then  the  approximate  variation  of  conductivity 
through  the  nozzle  is  obtained  by  simply  combining  the  equilibrium  variation  up¬ 
stream  of  the  freeze  point  with  a  straight  line  (on  the  log5  -  log  r  plane)  cor- 

P 

responding  to  the  frozen  portion  of  the  flow.  If  the  electron-neutral  collision 
cross-section  is  temperature  insensitive,  the  slope  of  these  lines  (represent¬ 
ing  lines  of  constant  electron  mole  fraction)  is  (Y  -  l)/V].  For  the  case  of 

Y  -  1.2,  the  slope  is  that  of  the  hypotenuse  of  the  triangle  depicted  in  Figs,  4* 
5  and  6.  This  general  procedure  is  described  in  Section  HI  of  the  text. 

Selection  o*f  Realistic  Parameters 


A  discussion  of  the  choice  of  the  exponents  and  is  in  order  since 
the  conclusions  depend  quite  markedly  on  the  choice  ultimately  made.  It  is  re¬ 
called  that  these  parameters  describe  the  dependence  of  the  recombination  coeffi¬ 
cient  on  pressure  and  temperature  and,  as  such,  reflect  differences  in  the 
microscopic  mechanism  of  ion  recombination.  Several  self-consistent  sets  are 
included  in  Table  VI,  together  with  remarks  on  the  applicability  of  each. 

For  a  fixed  temperature  level  these  mechanisms  are  such  that  each  may  be¬ 
come  dominant  at  an  appropriate  pressure  level.  This  accounts  in  part  for  the 
shape  of  such  recombination  coefficient  curves  as  shown  in  Fig.  11.  The  presence 
of  pressure  independent  mechanisms  (such  as  the  last  two  given  in  Table  VII)  would 
have  the  effect  of  introducing  pressure  independent  lower  bounds  on  not  shown 
in  Fig.  11.  It  should  also  be  remarked  that  the  pressure  regimes  at  which  these 
••exchanges”  of  mechanism  occur  themselves  shift  considerably  with  changes  in  the 
temperature  level.  This  fact  appears  to  be  extremely  important  in  rocket  motor 
applications  since  statements,  such  as  **the  Thompson  mechanism  should  not  apply 
at  pressures  over  one  atmosphere**,  often  leave  unstated  the  temperature  level  at 
which  this  conclusion  is  justified.  At  rocket  motor  temperatures  the  Thompson 
mechanism  may  actually  be  applicable  to  100  atmospheres,  but  this  is  an  unanswered 
question  for  experiment  to  decide.  It  is  also  true  that  while  the  pressure  range 
of  interest  in  rocket  motor  applications  may  appear  to  pass  through  one  of  the 
transition  regions  on  an  isothermal  plot  of  vs  pressure,  this  may  not  be  the 
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case  in  the  actual  expansion  because  of  the  corresponding  temperature  changes 
(c.f.  trajectory  shown  on  Fig.  11  for  Y  =  1.2)  in  isentropic  nozzle  flow. 


TABLE  VII 

PRESSURE  AND  TEMPERATURE  DEPENDENCE  OF  THE  ION  RECOMBINATION  COEFFICIENT»^>*^>*^ 


Mechanism 

Theory 

U) 

p 

Remarks 

Three-Body 

J.J.,  G. P.  Thompson 
(1924) 

+l-»  0 

Thermal  velocity  de¬ 
termines  rate  of 
successful  encounters 

Three  Body- 

Mobility 

Controlled 

P.  Langevin 
(1903) 

-1 

Thermal  velocity  ne¬ 
glected  compared  with 
that  due  to  field  in¬ 
duced  drift.  May  not 
apply  to  electron-ion 
recombination 

Two -Body} 
Dissociative 

D.R.  Bates 
(1950) 

0 

May  apply  to  monatom-- 
ic  gases  due  to  ion 
clustering 

Radiative 

Strueckelberg  and 
Morse  (1930)  etal. 

0 

a  <  10*'^  for  ion- 
r 

electron  recombina¬ 
tion  at  ambient  tem¬ 
perature 

'Spontaneous ' 
(dielectronic) 

0 

a  ^  10*^  to  10*^* 
r 

at  ambient  tempera¬ 
ture  for  ion-electron 
recombination 

As  described  in  the  text,  calculations  have  been  carried  out  for  both 
the  Thompson  and  Langevin  mechanisms  and  it  is  found  that  for  comparable  values 
of  the  recombination  parameter  C^,  the  Thompson  mechanism  predicts  much  larger 
deviations  from  equilibrium  (i.e. >  earlier  freeze).  This  is  not  surprising,  par¬ 
ticularly  in  view  of  the  important  difference  in  pressure  dependence  (see  Table 
VII).  Interestingly  enough,  the  Langevin  values  of  and  are  just  those  re¬ 
quired  to  annihilate  the  explicit  pressure  and  temperature  dependence  of  the 
first  term  on  the  right  hand  side  of  Eq.  (B-8)j  i .e. .  for  this  case 
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i  =  ■  c.  ='’[5^  -  5^3  -  i  3  ?■’  4  (B-W) 

However»  this  is  more  of  a  curiousity  than  an  essential  simplification. 

Accuracy  of  the  Sudden  Freeze  Approximation?  Improved  Approximations 

While  the  "sudden  freeze"  technique  lends  itself  to  the  presentation  of 
the  effects  of  ionization  nonequilibrium  in  a  compact  form  suitable  for  rapid 
graphical  calculation  (see  Figs.  4«5«6),  where  accuracy  is  importantj  improved 
solutions  will  probably  be  necessary.  For»  in  reality,  the  transition  from  near 
equilibrium  to  highly  nonequilibrium  flow  will  take  place  over  a  range  of  a  Mach 
number  and  strict  constancy  of  the  electron  mole  fraction  will  never  be  achieved. 

It  is  then  of  interest  to  investigate  the  nature  of  the  exact  solutions  to  Eq. 
(B-8),  both  to  determine  the  absolute  error  in  a  "sudden  freeze"  approximation 
as  well  as  to  develop  more  accurate  approximate  procedures. 

The  nature  of  the  exact  solution  is  readily  appreciated  by  reference  to 
the  behavior  of  the  integral  curves  of  Eq.  (B-8),  sketched  in  Fig.  13.  One 
notices  that  over  a  wide  Mach  number  range  only  small  departures  from  equilibrium 
are  "permitted"  but,  ultimately,  the  integral  curves  turn  away  from  the  equilibrium 
line.  Having  accomplished  this,  they  assume  an  analytic  form  readily  obtained 
from  Eq.  (B-8)  by  setting  <3*  »  c3|q,  i.e. 

^  =  -  c  (5)  »  (T)  ^  (3)-’  (5)» 

dx  r  T  ox 

But  this  equation  can  be  integrated  to  obtain 

U)  +1  U).-l 

•  P  f  ^ 

- - - -  dJ  (B-27) 

u 

o'  0  o 

where  the  subscript  o  is  used  to  denote  any  point  in  the  nozzle.  Since  Mach  num¬ 
ber  is  a  more  convenient  independent  variable,  it  is  useful  to  introduce  the 
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function 


t 


U)  +1 

P  ^ 


t  dx  oM 


in  terms  of  which  Eq.  (B-27)  may  be  written  (for  M  >  l) 


(B-28) 


h(m  )]  (B-a9) 

Wt 


These  integral  curves  of  Eq.  (B-26)»  which  make  no  reference  to  5  (and,  hence, 

eq 

are  independent  of  ©/T^),  are  sketched  in  Fig.  14.  A  study  of  this  figure,  to¬ 
gether  with  the  relation  between  Eqs.  (B-8)  and  (B-26),  now  suggests  an  alternate 
approximation  method  expected  to  be  more  accurate  than  the  "sudden  freeze"  tech¬ 
nique  described  earlier.  This  improved  approximation  method,  which  may  be  called 
a  "maximum  lower  bound"  technique,  also  patches  two  solutions  together  at  an  appro¬ 
priate  match  point  and  is  described  below. ^ 

The  relation  between  the  slope  fields  of  Eqs.  (B-8)  and  (B-26)  is  such  that 
the  result  of  following  the  equilibrium  curve  to  some  point  and  then  proceeding 
along  an  integral  curve  of  Eq.  (B-26)  will  always  fall  beneath  an  exact  solution 
of  Eq.  (b-8).  Clearly  then  the  maximum  lower  bound  will  everywhere  be  obtained  by 
following  the  equilibrium  curve  until  it  becomes  tangent  to  an  integral  curve  of 
Eq.  (B-26)  and  then  following  this  latter  integral  curve  beyond  the  tangency  point. 
At  this  tangency  or  "match"  point  (subscript  m)  both  5  and  do/dM  are  continuous, 
so  that 


_  _  "K+O 

u(p)  ^ 


(f)  "  0.,)- 


f/  ffal  .  i  ^  (Si 

Ll  dJ  1  f  dx 


(B-30) 


Note  that  this  is  identical  with  Eq.  (B-14)  (the  previous  "freeze  criterion"  with 
the  factor  ^  removed.  Riysically,  it  can  be  shown  that  m  represents  the  point 
on  the  equilibrium  curve  at  which  the  equilibrium  rate  of  change  of  electron  mole 

^Actually  H  should  be  written  h(Mj  Y,  '*'p)  • 

^This  improved  method  is  equivalent  to  that  suggested  by  Smith  in  Ref.  l6. 
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fraction  is  equal  to  that  produced  by  the  recombination  reaction  alone.  Fig.  15 
compares  the  solution  obtained  in  this  way  with  both  the  exact  solution  and  a 
"sudden  freeze"  solution.  While  the  latter  predicts  the  exit  conductivities  to 
within  better  than  a  factor  of  two,  it  is  much  less  successful  than  the  maximum 
lower  bound  method.  Since,  in  practice,  the  maximum  lower  bound  method  is  not 
much  more  difficult  to  apply,  H  is  to  be  recommended  for  future  nonequilibrium 
flow  studies  of  this  type.  Rapid  calculation  using  this  method  will  be  facili¬ 
tated  by  the  combined  use  of  Fig.  16  and  Table  VIII.  Fig.  16  shows  the  position 
of  the  match  points  m  as  determined  from  Eq.  (B-30)  for  two  distinct  recombination 
mechanisms.  Numerical  values  of  H(M)  are  given  in  Table  VIII,  to  be  used  in  con¬ 
junction  with  Eq.  (B-29). 
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TABLE  VIII 


THE  FUNCTION  H(M)  FOR  Y  = 

1.2  ISENTROPIC  EXPANSIONS  IN  A  O  = 

15°  HALF 

ANaE  CONE 

M 

H(M)^ 

T 

P 

3-Body* 

Dissoc,^^ 

1.0 

0.9535 

0.5645 

0.4461 

0.0000 

0.0000 

1.2 

0.9349 

0.0169 

0.0273 

1.4 

0.9144 

0.3417 

0.0482 

0.0861 

1.6 

0.8923 

0.2547 

0.0796 

0.1568 

1.8 

0.8691 

0.1856 

0.1063 

0.2305 

2.0 

0.8452 

0.1328 

0.1274 

0.3025 

2.2 

0.8209 

0.09363 

0.1433 

0.3707 

2.4 

0.7966 

0.06526 

0.1549 

0.4343 

2.6 

0.7724 

0.04512 

0. 1634 

0.4928 

2.8 

0.7487 

0.03102 

0.1694 

0.5464 

3.0 

0.7255 

0.02126 

0.1736 

0.5951 

3.2 

0.7029 

0.01455 

O.O29957’’ 

0. 1767 

0.6395 

3.4 

0.6810 

0.1788 

0.6797 

3.6 

0.6600 

O.O26826 

0.1803 

0.7162 

3.8 

0.6397 

0. 0,4692 

0.1813 

0.7494 

4.0 

0.6202 

0.0,3237 

0.1821 

0.7795 

4.2 

0.6015 

0.0,2243 

0.1826 

O.8O69 

4.4 

0.5836 

0.0,1561 

0.1830 

0.8318 

4.6 

0.5665 

0.0,1092 

0.1833 

0.8545 

4.8 

0.5501 

O.O37687 

0.1835 

0.8752 

5.0 

0.5345 

0.0j5440 

0.1836 

0.8942 

absolute  values  subject  to  +2  digit  maximum  error  in  last  place  tabulated 

*  Thompson  mechanism}  =+l,  =  -  7/2 

^^Dissociative  mechanism;  w  =  0,  =  -  3/2 

P  T 

^*The  subscript  appearing  on  the  zero  to  the  right  of  the  decimal  point  indicates 
the  number  of  zeros  appearing  before  the  first  non-zero  digit»  i .e.  0, 0^9957  r 

0.009957. 
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APPENDIX  C 

DECAY  OF  ELECTRICAL  CONDUCTIVITY  DOWNSTREAM  Cff  THE  NOZZLE 

In  the  region  beyond  the  nozzle  the  rate  of  change  of  stream  tube  area 
becomes  negligible,  so  that  all  gasdynamic  properties  of  the  working  fluid  tend 
to  remain  constant.^  However,  if  a  departure  from  ionization  equilibrium  has 
occurred  within  the  nozzle,  the  electron  recombination  reaction  (b-2)  persists 
and  it  is  of  interest  to  find  the  spatial  rate  of  approach  to  ionization  equili¬ 
brium  for  the  seed  material.  If  the  distance  required  to  return  to  equilibrium 
is  much  larger  than  the  contemplated  streamwise  extent,  L,  of  the  magnetic  field 
(see  Fig.  1)  then  no  significant  decay  of  the  electrical  conductivity,  CJ,  would 
occur  within  the  magnetic  interaction  region. 

In  a  downstream  region  of  constant  stream  tube  area,  Eq.  (B-8)  simplifies 
considerably  since 


p  =  constant  =  p 

T  =  constant  =  T 

0  „  =  constant  =  <3 
eq 

•A  ^  « 

u  =  constant  =  u 


ex 

ex 

eq,ex 

ex 


(c-l) 


As  a  result,  the  normalized  electrical  conductivity  o/d  will  decay  in  accord  with 

c 

the  equation 


dx 


U)  +1 

P 


^^ex^‘ 


[<5)'  -  (5.,..,)'] 


(C-2) 


In  this  region  it  is  therefore  convenient  to  introduce  two  new  variables,  the  first 
of  these  being  the  nonequilibrium  enhancement  ratio,  E,  defined  by 


S  = 


d/d 


ex,eq 


(C-3) 


For  a  specified  initial  value,  Z  will  depend  only  on  the  stretched  streamwise 
coordinate 


^apart  from  effects  such  as  radial  diffusion  and  turbulent  mixing  at  the  jet 
boundary 
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*  ^  (=ex'-’  '’.x.e.,  ''“Ax)  '-^x) 

Eq,  {C-2)  reveals  that  £(2)  must  satisfy  the  separable  first  order  equation 


E'  =  -  (Z*  -  1) 


(C-5) 


subject  to  the  initial  condition  Z(o)  =  Z^.  The  corresponding  solution  is  found 
to  be 

(Z  +  1)  +  (Z^  -  1)  exp  (-  2z) 

S(2)  =-^ - 2 -  (C-6) 

(Z^  +  1)  -  (Z^  -  1)  exp  (-  22) 

Thus  the  decay  of  nonequilibrium  electrical  conductivity  downstream  of  the  nozzle 
is  readily  calculable  in  terms  of  the  previous  nozzle  flow  parameters.  Results 
are  shown  in  Fig.  17.  Perhaps  of  greater  interest  in  the  present  study  is  the 
average  value  ^Z^  of  the  enhancement  ratio  Z  on  the  interval  from  0  to  z,  since 
this  average  would  be  more  representative  of  the  overall  interaction  region  than 
either  the  final  or  initial  values.  Fortunately,  a  simple  expression  for  the  aver¬ 
age  value  Z  is  obtainable  from  Eq.  (C-6),  viz. 

=  1  +  ^  In  +  1)  -  -  1)  exp  (-  2z)]  (C-7) 

This  relation  is  shown  graphically  in  Fig.  18  and  can  be  used  to  assess  the  impor¬ 
tance  of  electron  recombination  downstream  of  the  nozzle  in  decreasing  the  (spa¬ 
tial)  average  value  of  the  nonequilibrium  electrical  conductivity.  To  facilitate 
combining  these  results  with  those  of  the  nozzle  expansion,  analysis  of  Appendix 
B  plots  of  the  following  universal  function  (c.f.  Fig.  19)  are  useful 

+1 

G(p^/p  ;  0/T^,  Y.  <*>7*  ‘*'p)  =  P  ^  T  G'’  X  5^  (C-8) 

With  the  help  of  this  function,  the  downstream  distance  (measured  in  exit  diameters) 
is  readily  related  to  the  stretched  coordinate  z  as  follows 


z  =  •  G  •  (x/d^^)  (C-9) 

Here,  the  recombination  parameter  has  already  been  introduced  and  evaluated  for 
the  nozzle  flow  analysis. 
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To  illustrate  the  use  of  these  graphs^  consider  the  case  of  a  Y  =  1.2  isen- 

tropic  expansion  through  a  pressure  ratio  p  /p  =  68.  If  the  equilibrium  ioniza- 

c 

tion  parameter  0/T  is  15  and  the  recombination  parameter  C  has  been  estimated 

C  X 

to  be  10®  we  ask; 

(a)  What  electrical  conductivity  is  likely  to  prevail  10  exit  diameters 
downstream  of  the  nozzle  exit  section? 

(b)  What  is  the  corresponding  average  value  of  the  electrical  conductivity 
over  the  first  10  exit  diameters? 

For  simplicity  we  will  make  use  of  the  "sudden  freeze"  analysis*  assuming  that 
the  Thompson  three-body  recombination  mechanism  prevails.  Then  Fig.  6  is  appro¬ 
priate  and  we  find  that  5.^  =  0.039  and  5  „  =  .0223;  consequently  S  = 

©X  6X  j  ©Cj  O 

=  17*  Fig*  19  reveals  that  z  =0.0,17  C  •  (x/d  ),  so  that  10  dia- 

meters  downstream  (x/d^^  =  10)  we  have  z  =  (O.Ojl?) (lO®) (lO)  =  1.7  x  10"*.  The 

value  of  Z(z)/^  corresponding  to  z  =  1.7  x  10"*  and  =  17  is  found,  from  Fig. 

17,  to  be  about  0.78.  The  average  value  is  estimated  from  Fig.  18  to  be 

0.88.  Therefore,  the  corresponding  values  of  £  and  ^£^  are,  respectively, 

(0.78)(17)  =  13  and  0.88(17)  =  15.  From  the  definitions  of  £  and  <5:)>  and  the  fact 

that  5  ^  _  =  O.O223,  we  then  conclude  that  5  =  5  ^  =  0.0223(13)  =  .031  and 

ex^eq  ex^eq 

<&>  =3  ^  \ O  =  (0.0223)(15)  =  .034.  Thus,  the  electrical  conductivity  10 

ex^eq 

diameters  downstream  of  the  nozzle  exit  section  is  likely  to  be  about  3.1  percent 
of  the  electrical  conductivity  in  the  chamber.  The  average  conductivity  over  the 
first  10  diameter  length  is  approximately  3.4  percent  of  the  chamber  conductivity. 
By  way  of  comparison,  if  equilibrium  had  been  maintained  everywhere,  the  electri¬ 
cal  conductivity  at  the  10  diameter  location  would  have  been  only  0.23  percent  of 
the  electrical  conductivity  in  the  chamber.  As  illustrated  by  this  example,  one 
cannot  expect  good  precision  in  the  use  of  Figs.  17  and  18.  Where  this  is  impor¬ 
tant,  the  complexity  of  Eqs.  (C-6)  and  (C-7)  is  far  from  prohibitive  and  their 
direct  use  is  recommended.  In  general,  the  curves  presented  will  then  be  useful 
in  revealing  what  regimes  of  operation  will  be  encountered,  and  i4)ether  more  pre¬ 
cise  calculations  are  necessary. 
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FIG.  2  DEPENDENCE  OF  ATTAINABLE  TRANSVERSE  THRUST  ON  THE  MAGNITUDE  OF  THE 
MAGNETOGASDYNAMIC  INTERACTION  PARAMETER  Q 
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FIG.  4  THE  EFFECT  OF  ELECTRON-ION  RECOMBINATION  LAG  ON  EXIT  VALUES  OF  THE 
ELECTRICAL  CONDUCTIVITY  FOR  ISENTROPIC  NOZZLE  EXPANSIONS 

(Y  -  1.2|  Langevln  Recombination) 
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PRESSURE  RATIO,  p^/p 


FIG.  5  THE  EFFECT  OF  ELECTRON-ION  RECOMBINATION  LAG  ON  EXIT  VALUES  OF  THE 
ELECTRICAL  CONDUCTIVITY  FOR  ISENTROPIC  NOZZLE  EXPANSIONS 

(Y  =  1.2|  Dissociative  Recombination) 
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PRESSURE  RATIO,  pc/p 


FIG.  6  THE  EFFECT  OF  ELECTRON-ION  RECOMBINATION  LAG  ON  EXIT  VALUES  OF  THE 
ELECTRICAL  CONDUCTIVITY  FOR  ISENTROPIC  NOZZLE  EXPANSIONS 

(y  =  1.2|  Thompson  Recombination) 


53 


54 


FIG.  7  NOZZLE  SHAPE  TREATED  IN  RECOIBINATIOM  STUDY 
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FIG.  9  PRESSURE  DEPENDENCE  OF  EQUILIBRIUM  EXPANSION  FACTORj  (Y  =  1.2) 
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FIG.  10  PRESSURE  DEPENDENCE  OF  EQUILIBRIUM  EXPANSION  FACTORi  (Y  =  1.3) 
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FIG.  11  PRESSURE  DEPEKOENCE  OF  THE  ION  RECOMBINATION  COEFFICIENT  FOR  AIR 


RECOMBINATION  PARAMETER 
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FIG.  12  LOCATION  OF  THE  MACH  NUMBER  BEYOND  WHICH  THE  ELECTRON  MOLE  FRACTION 

IS  EXPECTED  TO  REMAIN  UNCHANGED 

fY  =  1.2|  Thomoson  Pocomblnatlon) 
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FIG.  13  EXACT  INTEGRAL  CURVES  OF  THE  DIFFERENTIAL  BQUAnON  GOVERNING 
CONDUCTIVITY  CHANGES  IN  THE  fCZZLE 
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FIG.  U  IMPROVED  APPROXIMATE  SOLUTIONSi  MAXIMUM  LOMER  BOUl®  TECHNIQUE 


FIG.  15  COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTIONS  FOR  A  PARTICULAR  CASE 
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FIG.  17  DECAY  OF  ELECTRICAL  CONDUCTIVITY  DOWNSTREAM  OF  THE  NOZZLE 


64 


62-52 


FIG.  18  DECAY  OF  SPATIALLY  AVERAGED  ELECTRICAL  CONDUCTIVITY  DOWNSTREAM  OF  THE  NOZZLE 
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